Ekpyrotic and Cyclic Cosmology 



Jean-Luc Lehners 



Princeton Center for Theoretical Science, Jadwin Hall, Princeton University, 

Princeton NJ 08544, USA 



Ekpyrotic and cyclic cosmologies provide theories of the very early and of the 
very late universe. In these models, the big bang is described as a collision of branes 
- and thus the big bang is not the beginning of time. Before the big bang, there 
is an ekpyrotic phase with equation of state w = ^ 1 (where V is the average 
pressure and p the average energy density) during which the universe slowly con- 
tracts. This phase resolves the standard cosmological puzzles and generates a nearly 
scale-invariant spectrum of cosmological perturbations containing a significant non- 
gaussian component. At the same time it produces small-amplitude gravitational 
waves with a blue spectrum. The dark energy dominating the present-day cosmo- 
logical evolution is reinterpreted as a small attractive force between our brane and 
a parallel one. This force eventually induces a new ekpyrotic phase and a new brane 
collision, leading to the idea of a cyclic universe. This review discusses the detailed 
properties of these models, their embedding in M-theory and their viability, with 
an emphasis on open issues and observational signatures. 
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1 Introduction 

The evolution of our universe is very well understood from the time of big bang 
nucleosynthesis until the present. However, if we go beyond these time fron- 
tiers, almost nothing is known with a great amount of certainty. We know that 
the early universe must have been in a very special state (very homogeneous 
and flat, but with tiny fluctuations in curvature), but we do not know why. 
The theory of inflation [41,69,4] provides a possible explanation by means of 
a period of rapid expansion preceding nucleosynthesis, but it has not yet been 
possible to test inflation conclusively, which is why it is important to keep an 
open mind. Also, we know that dark energy has come to dominate the energy 
density of the universe [86,85] and that it will determine the expansion of the 
universe in the near future. If dark energy turns out to be a cosmological con- 
stant, which is the theoretically simplest way of modeling it, then the universe 
will continue to expand at an accelerating pace and become increasingly dilute 
and cold, making life very difficult, if not impossible. However, it should be 
borne in mind that we have no fully convincing explanation of dark energy at 
present. 

This review will be concerned with a set of ideas, strongly inspired by string 
theory, that suggests alternative solutions to the early universe puzzles men- 
tioned above, and considers an alternative fate for the future of our universe. 
What makes these ideas exciting is that, on the one hand, they provide a the- 
oretical playground for applying string theory to cosmology, and, on the other 
hand, certain observational signatures predicted from these models are in a 
range that will be tested by near-future satellite experiments. This provides a 
certain timeliness to a subject concerned with eras far removed from everyday 
experience. 
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orbifoid dimension 

Fig. 1. The braneworld picture of our universe. Think of a sandwich: the fihing is the 
5-dimensinonal bulk spacetime, which is bounded by the two pieces of bread a.k.a. 
the 4-dimensional boundary branes. There is no space "outside" of the sandwich, 
but the branes can be infinite in all directions perpendicular to the line segment. 
In the M-theory embedding, there are 6 additional internal dimensions at each 
5-dimensional spacetime point. 

Ekpyrotic ^ and cyclic cosmology are based on the braneworld picture of the 
universe, in which spacetime is effectively 5-dimensional, but with one di- 
mension not extending indefinitely, but being a line segment ^ , see Fig. 1. The 
endpoints of this line segment (orbifoid) are two (3-|-l)-dimensional boundary 
branes. All matter and forces, except for gravity, are localized on the branes, 
while gravity can propagate in the whole spacetime. Our universe, as we see it, 
is identified with one of the boundary branes and, as long as the branes are far 
apart, can interact with the other brane only via gravity. The ekpyrotic model 

^ The name ekpyrosis can be translated as all-engulfing cosmic fire. In Stoic phi- 
losophy, it represents the contractive phase of eternally-recurring destruction and 
re-creation [1]. 

^ This setting will be discussed in detail in section 6, along with its motivation from 
and embedding in heterotic M-theory. 
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assumes that there is an attractive force between the two branes, which causes 
the branes to approach each other. This ekpyrotic phase has the rather non- 
intuitive property that it flattens the branes to a very high degree. Eventually 
the two branes coUide and move through each other (since there is no space 
"outside" of the boundary branes, it makes no difference whether we say that 
the branes bounce off each other or move through each other) . It is this colli- 
sion that, from the point of view of someone living on one of the branes, looks 
like the big bang. The collision is slightly inelastic and produces matter and 
radiation on the branes, where the standard cosmological evolution now takes 
place. However, due to quantum fluctuations, the branes are slightly rippled 
and do not collide everywhere at exactly the same time. In some places, the 
branes collide slightly earlier, which means that the universe has a little bit 
more time to expand and cool. In other places, the collision takes place slightly 
later, and those regions remain a httle hotter. This provides a heuristic picture 
of the way temperature fluctuations are naturally produced within the model. 
Shortly after the branc collision, the distance between the boundary branes 
gets almost stabilized, but the branes start attracting each other again very 
slightly. This very slight attraction acts as quintessence, and is identified with 
the dark energy observed in the universe. After a long time and as the branes 
become closer again, they start attracting each other more strongly so that 
we get a new ekpyrotic phase and eventually a new brane collision with the 
creation of new matter. In this way, a cyclic model of the universe emerges. 

The preceding description of course only provided a very rough outline of 
the main ideas. All of this will be made precise in the following. However, 
at this stage it is already clear that conceptually the cyclic universe differs 
substantially from the standard big bang picture. For one, the big bang is seen 
as a physical event and not a mysterious moment of creation. As such, it does 
not represent the beginning of time (note that if quantum gravity is unitary, 
this point of view seems unavoidable). Thus there was plenty of time before 
the big bang for the universe to be in causal contact over large regions, and 
in this way the horizon problem is automatically solved. The long timescales 
existing before the big bang also imply that other cosmological puzzles don't 
have to be solved in the extremely short time interval between the big bang 
and nucleosynthesis. What's more, properties of the very early universe (such 
as cosmological perturbations) and properties of the late universe (such as 
dark energy) are described by the same ingredients, namely the motion of 
branes. 

In the description above, the higher-dimensional interpretation was empha- 
sized. This higher-dimensional viewpoint is useful for understanding the ori- 
gin of the ideas, and provides a concrete realization of them. However, it is 
important to realize that many of the processes mentioned above can also be 
discussed purely in a 4-dimensional effective theory. Some processes, such as 
the ekpyrotic period of slow contraction before the big bang, are quite gen- 
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eral and could have a different origin than the brane motion just described. 
Therefore, the four subsequent sections of this review, which describe the fun- 
damentals of ekpyrotic and cyclic cosmology, are phrased mostly in terms of 
a 4-dimensional effective theory. Section 2 describes the ekpyrotic phase and 
the way it resolves the standard cosmological puzzles. Section 3 describes the 
approach to the big crunch in more detail. The cosmological fluctuations gen- 
erated by the model are discussed in section 4. This includes a treatment of 
scalar perturbations with their second-order non-gaussian corrections as well 
as tensor perturbations. The following section is concerned with the details of 
the cyclic universe, such as its relationship to the second law of thermodynam- 
ics and its global structure. Section 6 explains the origin of the 4-dimensional 
effective theory used in the preceding sections and its relation to string and 
M-theory. Finally, section 7 provides a brief conclusion and an outlook to 
possible future developments. 

It should be noted from the outset that ekpyrotic and cyclic cosmology are 
theories that arc still in full development. Certain aspects of the theories have 
gained or lost in importance over time. Here are two examples: in the original 
ekpyrotic proposal, the brane collision was between a boundary and a bulk 
brane, the second boundary brane being a spectator. It was later realized that 
the model could be simplified by considering a collision of the boundary branes. 
It was also thought that the dark energy phase was crucial for smoothing and 
isotropizing the universe. Again, it was later realized that the ekpyrotic phase 
can achieve this on its own. In this review the subject will be discussed only 
according to its present understanding. Moreover, every effort will be made to 
highlight assumptions that are in need of a precise calculation, or that are due 
for example to limitations in our understanding of string theory. This should 
hopefully provide the reader with opportunities to contribute to the subject. 



2 The Ekpyrotic Phase 



The ekpyrotic phase is the central insight that shows how a contracting phase 
preceding the big bang can solve the standard cosmological puzzles [48,35]. 
This is a surprising statement since one would naively think that a contracting, 
gravitating system would lead to large curvatures near a singularity. And we 
know from the near-flatness of our current universe that shortly after the big 
bang the universe must have been extremely fiat. This is usually referred to 
as the flatness problem. Let us briefly quantify this problem; at the same time 
this will serve to set up notation. Consider a Priedmann-Robertson- Walker 
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(FRW) metric^ 



ds^ = -dt^ + a{tf + r^dnl (1) 

\ 1 — nr-^ J 

where a{t) denotes the scale factor of the universe and n = —1,0,1 for an 
open, flat or closed universe respectively. If we consider the universe to be 
filled with a perfect fluid (which is a good approximation for many types of 
matter), with energy- momentum tensor 

T!^ = dtag[-p,V,V,V], (2) 

where p denotes the energy density and V the pressure, then the Einstein 
equations reduce to the so-called Friedmann equations 



^^4"-^ (3) 
^ = -i(p + 3P), (4) 

where the Hubble parameter is defined by if = - and a dot denotes a derivative 
w.r.t. time t. The Bianchi identity on the Einstein tensor VGfj_i, = requires 
for consistency that the energy-momentum tensor be covariantly conserved, 
i.e. V^T^i, = 0, which leads to the continuity equation 

p + 3^{p + V)^0. (5) 



ill Cl^UClLUJll UL iSLClLC UU = — 

can be integrated to give 



For a fluid with a constant equation of state w = ^, the continuity equation 



poca-'(^+"'). (6) 
We will use this scaling relation on many occasions. Now define the quantity 
Q{t) = j^^^ where pcru — denotes the critical density. Then the first 
Friedmann equation (3) can be rewritten as 

and it expresses how close the universe is to flatness. At the present time, 
observations show that [52] 

|Q-l|o<10-^ (8) 
Extrapolating back in time, this means that at the Planck time 

|Q-l|p; {aH)l 



\n-l\o {aHfpi 

r-2 



(9) 



I will mostly use natural units h = c=l and SttG = Mp, = 1. 
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If we assume a radiation-dominated universe {w = 1/3). which is a good 
approximation for the present calculation, then the continuity and Friedmann 
equations lead to p oc a~* and a oc t^/^, so that 

|l^-l|pz<^10-'~10-«'. (10) 

Even though extrapolating all the way back to the Planck time is probably 
exaggerated, this simple estimate shows that the universe must have been 
extremely flat at early times. Clearly, this peculiar observation asks for an 
explanation. 

In order to understand how the ekpyrotic phase can resolve this problem, it is 
useful to briefly review the solution proposed in inflationary theory. Inflation 
postulates a period of rapid expansion immediately following the big bang, 
during which the equation of state ty = ^ of the universe is close to — 1 
[41,69,4]. Prom the acceleration equation 

- = -^(p + 3^)^-^(1 + 3^), (11) 
a b 6 

we can see that such an equation of state indeed leads to accelerated expansion. 
One way to model a matter component with the required equation of state 
is to have a scalar field with canonical kinetic energy and with a very flat 
potential V{(f)), i.e. we add the following terms to the Lagrangian 

^(^-lg^'^d,(l>d.(l>-V{(l>)y (12) 

Then a quick calculation shows that (in the absence of spatial gradients) the 
equation of state is given by 

which is close to —1 if the potential is sufficiently flat so that the fields are 
rolling very slowly, i.e. if ^ ^(0)- Iii fhe presence of different matter 
types, represented here by their energy densities p, the Friedmann equation 
(3) generalizes to 



The subscript m refers to non-relativistic matter and includes dark matter, 

r refers to radiation and a denotes the energy density of anisotropics in the 
curvature of the universe (the associated scaling with a will be calculated in the 
following section). In an expanding universe, as the scale factor a grows, matter 
components with a slower fall-off of their energy density come to dominate. 
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<t>beg 



Fig. 2. The potential during ekpyrosis is negative and steeply falling; it can be 
modeled by the exponential form V{(j)) = —Voe~'^^. 

Eventually, the inflaton, whose energy density is roughly constant, dominates 
the cosmic evolution and determines the (roughly constant) Hubble parameter 
while causing the scale factor to grow exponentially, 

a oc e"K (15) 

We can define the relative energy density in the curvature as = —n/{o?'H'^) 
and in the anisotropies as fi^r = cr^/(3a^if^). During inflation, these relative 
densities fall off quickly, and the universe is rendered exponentially flat; ac- 
cording to (10) the flatness puzzle is then resolved as long as the scale factor 
grows by at least 60 e-folds. 

Now we will show how the same problem can be solved by having a contracting 
phase before the standard expanding phase of the universe. The Friedmann 
equation relates the Hubble parameter to the total energy density in the uni- 
verse, which is the sum of kinetic and potential energy. Now suppose that, 
instead of a flat potential, the scalar has a very steep, negative potential, as 
shown in Fig. 2. As a concrete example, one can model the potential with a 
negative exponential 

r(0) = -Voe-'t (16) 

where Vq and c are constants. Let us give one brief motivation for a negative 
potential: foreshadowing a cyclic picture, we can ask if it is possible for the 
universe to revert from an expanding to a contracting phase in a non-singular 
way. Then at some point the Hubble parameter must go through zero, which 
can be achieved by having a negative potential in order to cancel the positive 
kinetic energy of matter. From (13) this automatically implies an equation 
of state w > 1, and will be seen to have very powerful consequences. In a 
contracting universe, the argument presented in the previous paragraph is 
reversed, and one would initially expect the anisotropy term (proportional to 
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a~^) to come to dominate the cosmic evolution. However, if there is a matter 
component with w > 1, then (6) imphes that this component will scale with an 
even larger negative power of a, and hence will come to dominate the cosmic 
evolution in a contracting universe in the same way as the inflaton comes to 
dominate in an expanding universe. 

In fact it is straightforward to generalize the treatment to having many scalars 
(pi with potentials Vi(0i). Then, in a flat FRW background, the equations of 
motion become 

0i + 3H^i + Vi,^, = (17) 

and 

where Vi^^^ = {dVi/dcpi) with no summation implied. A useful relation is 
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(18) 



H--\T.k'. (19) 

If all the fields have negative exponential potentials Vi{(pi) = —Vi e"^^'^^ and if 
Ci ^ 1 for all i, then the Einstein-scalar equations admit the scaling solution 



Thus, we have a very slowly contracting universe with (constant) equation of 
state 

to = |- - 1 » 1. (21) 

We are using a coordinate system in which the big crunch occurs at t = 0; 
in other words, the time coordinate is negative during the ekpyrotic phase. 
The steeply faUing scalar fields act as a very stiff fiuid, and, in fact, one can 
take the condition w ^ 1 (or equivalently p ^ |) to be the defining feature 
of ekpyrosis. The matter content docs not necessarily have to be composed of 
scalar fields, but it is easy to model the ekpyrotic phase that way, and scalar 
fields commonly appear in effective theories arising from higher dimensions. 
The main consequence is that the extra term in the Priedmann equation (14) 
with w ^ \ comes to dominate the cosmic evolution, and once more, the 
fractional energy densities oc a~^H~^ and fi^ oc a~^H~'^ quickly decay. 
Thus, neglecting quantum effects, the universe is left exponentially flat and 
isotropic as it approaches the big crunch. As we will see in section 4, the 
inclusion of quantum effects superposes small fluctuations on this classical 
background. 

What are the requirements on the ekpyrotic phase for the flatness problem to 
be solved? First, we must ensure that w > 1, which means that the potentials 
have to be steep enough. In the single-field case for example, from (21) we see 
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that wc need Ci > \/6- Also, from the scahng solution (20), we see that the 
scale factor a remains almost constant, while the Hubble parameter H cc 
during ekpyrosis. In order to solve the flatness problem, aH has to grow by at 
least 60 e- folds in magnitude, see (7) and (10), which means that 

\tbeg\ ^ e^'^l^endl, (22) 

where the subscripts beg and end refer to the beginning and the end of 
the ekpyrotic phase respectively. As will be discussed in section 4, we need 
tend ~ — lO^Mp;^ in order to obtain the observed amplitude of cosmological 
perturbations, so that 

\tbeg\ > IQ^^Mpl = IQ-^^S. (23) 

This is the minimum time the ekpyrotic phase has to last in order to solve 
the flatness problem, under the assumption that a flat universe undergoing 
a crunch/bang transition can reemerge as a flat universe afterwards. We will 
discuss this assumption repeatedly in what follows. 

In the cyclic picture of the universe, the potential V interpolates between the 
GUT scale and the dark energy scale. From the scaling solution (20), we have 
that V oc t^^, and this leads to 



\tbeg\ = Jirf^ltendl ~ VW^lO'Mpl = lO^S. (24) 
V ^beg 

In this case, the ekpyrotic phase lasts on the order of a billion years, and it 
easily resolves the flatness puzzle. 

From the higher- dimensional point of view, one of the scalar fields 0, is 
the radion, which determines the distance between the branes (the higher- 
dimensional picture will be discussed in detail in section 6). The potential then 
represents a conjectured attractive force between the end-of-the-world branes. 
By diluting inhomogeneities in the brane curvature, the ekpyrotic phase causes 
the branes to become very flat and parallel over large patches. Thus the ho- 
mogeneity puzzle, namely the question of how different parts of the sky can 
have approximately the same temperature despite the fact that they could not 
have been in causal contact since the big bang, is also solved by the ekpyrotic 
phase. 

Moreover, as will be discussed in sections 3 and 6, the brane scale factors 
remain finite even at the collision. Semi-classical calculations confirm the in- 
tuition that the collision should be slightly inelastic, so that matter is produced 
at a finite temperature [104]. And as long as this temperature is below the 
GUT scale (which might depend on the collision rapidity being low enough), 
topological defects will not form, and thus the monopole problem is avoided. In 
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this way all the standard cosmological puzzles can be resolved by the ekpyrotic 
phase. 



3 Towards the Crunch 

We have just seen how the simple assumption of a contracting phase with 
equation of state w > 1, preceding the usual expanding phase of the universe, 
manages to solve the flatness and homogeneity puzzles of standard big bang 
cosmology But now we will see that it does much more than that: it provides 
a way of avoiding Belinsky-Khalatnikov-Lifshitz (BKL) chaotic mixmaster be- 
havior [12], it prepares the universe in the most favorable state as it approaches 
the big crunch and, on top of that, it generates a nearly scale- invariant spec- 
trum of scalar perturbations that act as the seeds for the large-scale structure 
of the universe. 

3.1 Avoiding Chaos 

Any cosmological theory involving a contracting phase must address the issue 
of chaos: from the studies of BKL [12], it is known that if all matter compo- 
nents have an equation of state w < 1, then a contracting universe is unstable 
w.r.t small perturbations. The metric becomes highly anisotropic, and of Kas- 
ner form (see below). Typically all spatial dimensions except for one shrink 
away, but repeatedly the metric jumps from one Kasner form to another one, 
thus leading to BKL oscillations or "chaotic mixmaster" behavior [79] . In this 
case all predictability in going towards a big crunch seems to be lost. 

In order to better understand the mechanism of the ekpyrotic phase, and to 
address both the issues of chaos and cosmological perturbations, it is useful 
to look in more detail at the evolution of anisotropics. In synchronous gauge, 
we can write the metric as 



In this gauge, we can choose a time coordinate such that the big crunch 
occurs everywhere t — 0. Prom the work of BKL [12] we know that in a 
contracting universe spatial gradients quickly become irrelevant compared to 
time gradients, so that we can consider the simplified metric 



which is of Kasner form and where we require Y^iPi — 0. In other words. 



ds^ = -di^ + gij{x'')dx'dx^ . 



(25) 



ds^ = -df + a(t)2 ^e2'3^WdxMx\ 



(26) 



11 



the dynamics becomes ultralocal and the Einstein equations reduce to the 
Friedmann equations 



(27) 



(28) 



Then the growing-mode solution to the last equation is /3j oc a 
that the first Friedmann equation can be rewritten as 



^, which means 




(29) 



for a constant a parameterizing the anisotropic energy density at unit scale 
factor. Hence, as stated in the previous section, anisotropics blueshift propor- 
tionally to and would quickly dominate the cosmological evolution if it 
wasn't for the ekpyrotic scalar whose energy density grows faster, namely as 
From the scaling solution (20) it is easy to see that the time-dependent 
part of the Kasner exponents /3i scales as (— 1)^~^^ so that these exponents 
approach constant values during the ekpyrotic phase where P < \- This gives 
rise to the so-called cosmic no-hair theorem for ekpyrosis [36], according to 
which a contracting universe with anisotropy and inhomogeneous curvature 
converges to a homogeneous, fiat and isotropic universe if it contains energy 
with equation of state w > 1. From (14) and (22) we can see that the rela- 
tive importance of the anisotropy term in the Friedmann equation drops by a 
factor of at least 



during the ekpyrotic phase (note that the relative importance of the curvature 
term oc a'^H'^ drops off by a similar amount). Thus, at the end of the 
ekpyrotic phase, the size of pre-ekpyrotic anisotropics and spatial curvatures 
has become exponentially small! 

The ekpyrotic phase does not last right up to the big crunch. Indeed, on 
general grounds one expects the potential to turn off and become irrelevant 
as the big crunch is neared, see section 6. Then the universe enters a kinetic 
energy dominated phase, with equation of state w — 1. The equations of 
motion reduce to 




(30) 



2 



-H 



+ 3H(I) = 0. 



(31) 



They immediately lead to a oc e' 



and are solved by 





(32) 
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for some integration constants oq. 0o- As will be derived in section 4, the 
ekpyrotic phase lasts until about 10'^ Planck times before the big crunch; 
the subsequent kinetic phase is relevant until the quantum gravity regime is 
reached around t ~ — IMp/, and also again in the first moments after the big 
bang. Now of course one has to answer the question of whether this kinetic 
phase undoes what the ekpyrotic phase has achieved. From the solution above, 
one can see immediately that 0^ oc a"^, which means that the scalar field 
energy density scales in the same way as the anisotropies. Thus the relative 
importance of the anisotropic term remains constant during the kinetic phase. 
In a similar way [36], one can show that the curvature terms are harmless 
during this w = 1 phase. In the presence of p-forms, depending on their 
coupling to the scalar 0, dangerous modes can grow, but only as a power of 
time t and for not longer than about 7 e-folds. In the relevant case of heterotic 
M-theory for example [73], one such coupling is critical [36] implying that 
the approach to the big crunch could be chaotic if there was sufficient energy 
density in the p-form modes. 

Now it is clear how the ekpyrotic phase can avoid the uncertainties of a chaotic 
crunch: even when dangerous modes are present their energy density gets 
diluted to such an extent during ekpyrosis that they don't have time to come 
to dominate before less than a Planck time before the big crunch*^ [36]. Thus 
chaos is avoided by being delayed until we have entered the quantum gravity 
regime, where the equations of general relativity break down. At that point, 
one might heuristically expect the spread of the wavefunction of the universe 
to prevent chaotic behavior from resurfacing. This last speculation will be 
discussed further in section 6. 

3.2 The Milne Universe 

After the ekpyrotic phase, the universe enters a kinetic energy dominated 
phase as it approaches the big crunch. Thus, from a 4-dimensional point of 
view, the universe seems to be headed towards the worst kind of singularity. 
However, in approaching the crunch at t = 0, the universe really starts looking 
5-dimensional, and the higher-dimensional description reveals that the singu- 
larity is actually much milder. We will give a simplified description here; a 
more complete treatment is postponed to section 6. 

Consider gravity in 5 dimensions, but assume that one spatial dimension is a 
line segment, as suggested by Hofava-Witten theory (see section 6). This line 
segment can be described as a circle modded out by a Z2 reflection symmetry 
across an axis, see Fig. 3. Then we can write the metric as [47] 

^ What's more, some dangerous modes are further suppressed by the topology of 
the internal manifold, see [107]. 
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Fig. 3. A line segment represented as the orbifold 3^/1^2- The thick dots represent 
the orbifold fixed points. 

dsl = e-v^-^ds^ + e^v^^dy^, (33) 

where — ?/o < y < Ho denotes the orbifold coordinate and is the radion field 
parameterizing the size of the orbifold. The ^-dependent prefactor in front of 
the 4- dimensional metric ds4 ensures that after dimensional reduction we are 
left with the canonically normalized Lagrangian 

c = V^[lR-l{d<pY]. (34) 

We know that the kinetic phase described in equation (32) is a solution to this 
theory. This means that we can now plug the solution (32) into the original 
5-dimensional metric, to find (up to trivial rescalings) 



dsl = (-t)-2/3[_dt2 + i-tf/Mxl] + i-ty/'dy' (35) 
= -dT^ + T^y^ + dxl (36) 

after changing coordinates to T oc (— t)^/'^. This spacetime goes by the name 
of "compactified Milne mod Z2" XM3 [99]. It describes the two orbifold planes 
{i.e. the (3 + l)-dimensional boundaries of the line segment) approaching each 
other, colliding, and receding away from each other again, as shown in Fig. 4. 
Under the further change of coordinates u = Tcoshy, v = Tsinhy, the metric 
becomes 

ds^ = -dw^ + dv"^ + dxl, (37) 

which shows that the embedding spacetime is simply Minkowski space. Thus, 
except at the collision itself, this spacetime is fiat and therefore automatically 
a solution of any theory of gravity without a cosmological constant. This 
also means that higher-derivative corrections are small for a spacetime that 
is close to Milne. In this sense the compactified Milne space is the simplest 
model spacetime for a brane collision. Note that, from the 5-dimensional point 
of view, the big crunch singularity is just the momentary shrinking away of 
the orbifold dimension, while the other dimensions remain at finite size. This 
can best be seen by examining equation (33) together with the solution (32): 
even though the 4-dimensional scale factor a — > at the collision, the brane 
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V 



Fig. 4. The compactified Milne mod Z2 space describes the colhsion of two boundary 
branes, embedded in Minkowski space. 

scale factors are really given by e~'^l^a^ which is a constant. Thus the density 
of matter (which is stuck on the branes) and the temperature remain finite at 
the collision. 

The general idea here is that the ekpyrotic phase drives the universe to being 
exponentially flat and isotropic in a 4-dimensional sense, which means that 
the subsequent kinetic phase and big crunch/big bang transition can be well 
approximated in 5 dimensions by the compactified Milne spacetime. Of course 
one can only trust the spacetime description up to about a Planck time on 
either side of the collision. For a proper treatment, we would need a complete 
theory of quantum gravity. But in a sense, the ekpyrotic phase prepares the 
universe in the best possible state to enter the collision phase. 

There have been a number of studies in string/M-theory of the Milne spacetime 
near the collision: Turok et al. [104] have shown how the equations for winding 
membranes remain well behaved near the collision. And in [82] it was argued 
that from the string theory perspective strings evolve as sets of weakly coupled 
"string bits" smoothly across the singularity. However it is not clear to what 
extent these treatments encompass all the relevant degrees of freedom near the 
collision. Furthermore, it was argued in [104] that near the collision the theory 
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should be described in terms of an expansion in the string tension \jci rather 
than in such an expansion is not well understood at present. It is also known 
that perturbative string theory breaks down in these and similar settings as 
soon as interactions are taken into account [71,13,27]. But again it is not 
clear that perturbative string theory is really the right approach. The most 
promising development so far is the treatment of a toy model crunch/bang 
transition in the framework of the AdS/CFT correspondence by Craps et al. 
[103,26]. However, in this toy model, the bulk spacetime is 5-dimensional and of 
Kasner type close to the crunch. But most importantly, the radius of curvature 
of the AdS spacetime is smaller than the string length, in sharp contrast to 
ekpyrotic models, in which we rather expect the curvature to be small in 
approaching the crunch. Hence the results of this work cannot directly be 
applied to ekpyrotic models. Nevertheless, one might hope that the techniques 
used in [26] can be extended in the future to include the treatment of more 
realistic spacetimes. 

Of course this is a topic of crucial importance for ekpyrotic theory. In order to 
discuss cosmological perturbation theory, it is vital to know how to match the 
universe going into the big crunch to the one coming out of the big bang. In 
light of the ambiguities discussed above, two broadly different attitudes can be 
adopted (see [10,20,29,34,45,74,76,77.98,99,102] for many discussions of these 
"matching conditions" in the literature): 

In analogy with the treatment of reheating in inflation, Creminelli et al. [29] 
have advocated the model-independent assumption that the spacetime metric 
can simply go through the collision essentially unchanged. This assumption 
is based on the following argument: the ekpyrotic phase renders the universe 
isotropic up to exponentially small terms. Also, as we will see in the next 
section, quantum effects lead to small anisotropies in the metric. Thus, in 
synchronous gauge (g^o/i = Vo/j.) , the metric reads 

ds^ = -dt^ + a(i)2e2C(-'"e^'^-(^'^)da;Ma;^ (38) 

where we are allowing for scalar curvature perturbations ( and tensor per- 
turbations hij. Then locally the small, long- wavelength perturbations we're 
interested in can be gauged away by a local re-scaling of the spatial coordi- 
nates. In other words, locally these fluctuations are just integration constants. 
Thus the universe undergoes the same history at every point in space up to ex- 
ponentially small terms. If we then assume that the dynamics is not sensitive 
to these exponentially small terms (on super-horizon scales) then the metric 
will reemerge from the crunch/bang transition with the same integration con- 
stants, i.e. with the same perturbations. Note that this argument, if correct, 
applies not only to linear perturbations but to the full non-linear metric. 

The second possibility is that the dynamics of the crunch/bang transition is 
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important: for example it was shown in [98,99] that for free fields in a com- 
pactified Milne universe, one can simply analytically continue the evolution 
of the fields around the singularity at t = 0. Such a prescription corresponds 
to matching perturbations on surfaces of constant energy density (see also 
[49,105]) and it leads to a mixing of different perturbation modes, such as the 
dominant and subdominant components of the Newtonian potential and the 
curvature perturbation. As the energy density perturbation is exponentially 
small, this represents a concrete example of a situation in which the history 
of the universe does depend significantly on the exponentially small perturba- 
tion terms. It should be noted however that this treatment breaks down once 
interactions are taken into account, as noted above. It has also been suggested 
by McFadden et al. [78] that a mixing of modes can occur when the effective 
4-dimensional description breaks down near the collision. It should be clear 
that in this case as well, the evolution of the metric and its perturbations 
depend sensitively on the detailed dynamics of the bounce. 

We will adopt the "nothing happens" assumption of Creminelli et al. in this 
review {i.e. we will assume that the curvature perturbation is non-linearly 
conserved through the collision), but with the proviso that this assumption 
is subject to revision once the crunch/bang transition is understood in more 
detail. 



3.3 New Ekpyrotic Models 

The idea behind the "new ekpyrotic" models [22,23,30] is to build cosmological 
models in which the ekpyrotic contracting phase is joined with the standard 
expanding phase by a smooth, non-singular bounce that can be described 
entirely within a 4-dimensional effective field theory. Thus, the energy never 
reaches the Planck scale, and cosmological perturbations can be evolved un- 
ambiguously through the boimce phase. Of course, in light of the singularity 
theorems of Hawking and Penrose, it is clear that in order to achieve a smooth 
bounce, one of the assumptions going into these theorems has to give; in the 
models considered so far this has been the null energy condition (NEC), which 
is violated by means of a ghost condensate [6] . ^ In the presence of a fluid with 
energy density p and pressure V, the Friedmann equations can be combined 
to yield 

H^-lip + V). (39) 

The null energy condition requires that p + V ^ 0, or, in covariant form, 
Ti^i^n'^n^ ^ for any non-spacelike vector n'^. Thus, in order to have H > 



^ We should note though that it is not clear yet if such a ghost condensate can be 
realized in string theory [3] or in a UV complete theory in general [46] . 
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0, which is a necessary condition to revert from a contracting phase to an 
expanding one, the NEC must be violated. 



We will briefly review the main ideas of ghost condensation, in order to see 
how the NEC can be violated: we start by considering an effective Lagrangian 
of the form 

C^^M^P{X), (40) 

where 

^ = -iW)" (41) 

and P{X) is a function that remains to be specified. M and m are two mass 
scales that must be determined by the underlying microscopic theory, and, as 
we will see, that have to satisfy certain consistency conditions. Note that the 
theory admits the constant shift symmetry i— > constant. In a cosmological 
background, the resulting equation of motion reads 

A(a3p^0)=O, (42) 

which is automatically solved at an extremum Px = 0. We will assume that 
the extremum lies at = | say, which in turn can be solved by giving the 
scalar the time-dependent expectation value 

= -mH. (43) 

With this expression for 0, we are justified a posteriori to have omitted higher- 
derivative terms like in the Lagrangian, since these would vanish in any 
case. 



The ghost condensate by itself does not violate the null energy condition yet; 
in fact it has the same equation of state as a cosmological constant. To see 
this, note that the energy-momentum tensor is given by 

T^, = M'Pg^, + —P^xd^4)d,4>- (44) 

Thus the energy density is given by p = M'^{2P^xX — P), while the pressure 
is P = M^P. Hence, at the extremum, we have V — —p or w — —1. However, 
the fluctuations around the extremum can violate the NEC. We deflne the 
fluctuations tt via 

= -mH + TT. (45) 
Their dynamics can be expressed by the effective Lagrangian 

C (X 2XoP,xx{Xo)n\ (46) 

which shows that the kinetic term has the correct sign if the extremum is a 
minimum, as we will assume (in other words, there are no ghosts around the 
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P(X), 



P(X) 



P(X) 





Fig. 5. This figure shows the evolution of the kinetic function P{X) and of tlic poten- 
tial V{(j)) during the phases of ekpyrosis and bounce. Note that the field starts at ap- 
proximately zero potential energy and moving very slowly, so that p ~ — vr + V ^ 0. 
Then, during the ekpyrotic phase, the potential is negative, implying that tt < and 
hence H < 0. Subsequently, during the ghost condensate phase, the potential shoots 
back up to positive values, implying that vr > and thus H > 0, i.e. the universe 
starts reverting from contraction to expansion. Figure reproduced with permission 
from [22]. 

ghost condensate). If we then expand the expression for the energy density, 
p — M^{2P^xX — P), to linear order in tt, we obtain p ^ —KM^-k/m?. We can 
provoke such a fluctuation by adding a potential V, implying that the energy 
density and the pressure are now approximately given by 



P 



■n + V V ^-V, 



(47) 



where K = P,xx{Xo) > and P{Xq) = 0. Prom (39) we can see that this 
immediately implies 

and hence, since tt can take either sign, we now have the possibility of violating 
the NEC. 



Nevertheless, it is still difficult to achieve a bounce because NEC-violating 
modes scale with a positive power of the scale factor, see equation (6). Hence 
they tend to become subdominant in a contracting universe, and dominant in 
an expanding one - exactly the opposite of what we want here! This problem 
may be circumvented if it is the same field driving both the ekpyrotic phase and 
the bounce. For this to work though, we need the Lagrangian function -P(X) 
to be canonical (linear) during ekpyrosis [22] (although see [30] for a different 
approach), and then quadratic around the minimum; also, the potential V 
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has to become positive at the onset of the bounce phase in order to push 
the equation of state parameter w to NEC- violating values below —1. The 
complete evolution is sketched in Fig. 5. After the bounce is complete, it is 
assumed that the NEC- violating fluid decays by converting all its energy into 
matter and radiation, and that it thereby reheats the universe [22]. 

In these models, there are two sources of instability, which must be addressed. 
One is of the Jeans type and there is also a gradient instability, see [22,28,30] 
for details. Both of these instabilities are harmless if the bounce phase is 
rather fast, i.e. if it lasts only on the order of one Hubble time (as measured 
at the start of the bounce phase). In that case, there isn't enough time for the 
instabilities to grow. A further constraint that these models have to satisfy 
in order to obtain a successful bounce is that the value of the potential at 
the end of the ekpyrotic phase has to fall within a certain range, so that the 
ekpyrotic phase lasts long enough and so that the trajectory ends up near 
the ghost-condensate point at the end of ekpyrosis. As discussed in [23], this 
bound is 

m^e^^* < iKndl < , 49 

V 

where Mek — \ ^'^\^endl^heg\ IS the number of e-folds of the ekpyrotic phase; 
in other words, we must require that M is exponentially larger than m. Of 
course it would be satisfying to see how this large hierarchy in the effective 
theory could arise from an underlying, more fundamental, theory. 

But this is not our direct concern here. The new ekpyrotic models rather aim 
to represent a proof-in-principle that it is possible to have a cosmological his- 
tory that can be described entirely within the framework of a 4-dimensional 
effective theory, so that the evolution of all fields, and in particular of cosmo- 
logical perturbations, can be followed unambiguously through the bounce. We 
will return to the issue of cosmological perturbations during the bounce phase 
in the next section. 



4 Cosmological Perturbations 

The ekpyrotic phase not only manages to address the standard cosmological 
puzzles, but it also provides a new way of generating cosmological pertur- 
bations. Both scalar and tensor perturbations are generated, but, as we will 
see, they turn out to have very different properties. The colliding branes pic- 
ture provides a heuristic way of understanding the origin of the perturbations 
[48]: although the ekpyrotic phase renders the branes very flat and parallel 
classically, quantum fluctuations cause the branes to ripple slightly. Thus, the 
branes do not collide everywhere at exactly the same time, and therefore, in 
some places the big bang happens slightly earlier than in others. In this way. 
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some places have a little bit more time to expand and cool than others, and 
the result is a pattern of temperature fluctuations in the cosmic microwave 
background (CMB). We will now make this cartoon quantitatively precise 



4 ■ 1 Scalar perturbations 

The detailed shape of the angular anisotropy power spectrum of the CMB 
radiation can be described by fluctuation modes that enter the horizon in 
succession, i. e. in a timed, coordinated manner. This hints to having a period 
of cosmological evolution where the modes were in causal contact initially, and 
then went outside of the horizon, only to re-enter later in succession. There 
are, broadly speaking, two ways of producing this phenomenon: in inflation, 
the horizon is roughly constant and the modes get stretched rapidly beyond 
it, whereas in ekpyrosis, due to the increasing ultralocality in approaching the 
crunch, the horizon shrinks rapidly while the fluctuation modes themselves 
remain roughly constant in size. As we will see, for scalar, linear perturbations 
this leads to a kind of duality between the two theories. Tensor modes, on the 
other hand, are sensitive only to the evolution of the metric, and hence behave 
very differently in the two theories; this will be discussed in section 4.2. 

4.1.1 Single Field 

We will start with a simple, heuristic argument for why one might expect to 
obtain scale-invariant perturbations in the presence of a scalar fleld with a 
steep negative potential [61]. In fact, since the universe is contracting very 
slowly during the ekpyrotic phase, to a flrst approximation we can ignore 
gravity altogether. The validity of this approximation will be discussed shortly. 
So, for now consider a scalar fleld (j) in Minkowski spacetime, with action 



We consider the case where Vq is positive so the potential energy is formally 
unbounded below and the scalar fleld runs to —00 in a flnite time. Note that 
the actual value of Vq is not a physical parameter, since by shifting the fleld 
one can alter the value of Vq arbitrarily. Furthermore, in the ekpyrotic or cyclic 
models, the potential is expected to turn up towards zero at large negative 
(because this corresponds to the limit in heterotic M-theory in which the 
string coupling approaches zero and the potential disappears [93]), but this 
detail is irrelevant to the generation of perturbations on long wavelengths. 

We now argue that, as the background scalar fleld rolls down the exponential 



[48,49,39,18]. 




(50) 
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potential towards — oo, then, to leading order in h, its quantum fluctuations 
acquire a scale-invariant spectrum as the result of three features. First, the 
action (50) is classically scale- invariant. Second, by re-scaling — > 0/c and 
re-defining Vq, the constant c can be brought out in front of the action and 
absorbed into Planck's constant h, i.e. /i — > h/c^, in the expression iS/h gov- 
erning the quantum theory. Finally, it shall be important that has dimensions 
of mass in four spacetime dimensions. 

To see the classical scale-invariance, note that shifting the field 0^0 + ^, 
and re-scahng coordinates x^^ — > x^e'^^'^, just re-scales the action by e*^^ and 
hence is a symmetry of the space of solutions of the classical field equations. 
Now we consider a spatially homogeneous background solution corresponding 
to zero energy density in the scalar field. In a cyclic model (see section 5), this 
"zero energy" condition is a reasonable initial state to assume for analyzing 
perturbations because the phase in which the perturbations are generated is 
preceded by a very low energy density phase with an extended period of ac- 
celerated expansion, like that of today's universe, which drives the universe 
into a very low energy, homogeneous state. No new energy scale enters and 
the solution for the scalar field is then determined (up to a constant) by 
the scaling symmetry: 06 = (2/c) ln(— ^4^). Next we consider quantum fluc- 
tuations S(f) in this background. The classical equations are time-translation 
invariant, so a spatially homogeneous time-delay is an allowed perturbation, 
4> = (2/c) lii{—A{t + 6t)) ^ 6(j) (X t~^. On long wavelengths, for modes whose 
evolution is effectively frozen by causality, i.e. \kt\ <^ 1, we can expect the per- 
turbations to follow this behavior. Hence, the quantum variance in the scalar 
field is 



Restoring c via (p —>■ c(f) and h ^ c^h leaves the result unchanged. However, 
since S(f) has the same dimensions as in four spacetime dimensions, it follows 
that the constant of proportionality in (51) is dimensionless, and therefore that 
Sep must have a scale-invariant spectrum of spatial fluctuations. 

It is straightforward to check this in detail. Setting (p = 4>b{t) + 50(t,x), to 
linear order in 50 the field equation reads 



Using the zero energy condition for the classical background, we obtain V^^^ = 



c^V = -c^(i)b /2 = -2/^2. Next we set (50(t,x) = Ek{a]^Xk{t)e^^"' + h.c.) 



with ajc the annihilation operator and Xk if) the normalized positive frequency 
modes. The mode functions Xk obey 



(502) (X hr''. 



(51) 



(52) 



^Xk = T^Xifc - k'^Xk, 



(53) 
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and the incoming Minkowski vacuum state corresponds to [14] 



(64) 



For large \kt\, this solution tends to the usual Minkowski positive- frequency 
mode. But as \kt\ tends to zero, each mode enters the growing time-delay 
solution described above, with Xk oc t"^. Computing the variance of the quan- 
tum fluctuation and subtracting the usual Minkowski spacetime divergence, 
we obtain 

where the integral is taken over k modes which have "frozen in" to follow 
the time-delay mode. In agreement with the general argument above, we have 
obtained a scale- invariant spectrum of growing scalar field perturbations. To 
recap, classical scalc-invariance determines the t dependence of the perturba- 
tions, and dimensional analysis then gives a scale-invariant spectrum in /c, in 
three space dimensions. 

However, it is important to realize that what is measured in the CMB is di- 
rectly related to the curvature perturbation of the spacetime geometry. Hence 
it is imperative that we include gravity in the analysis, even though a priori its 
effects appear to be small corrections. As we will see, the inclusion of gravity 
entails some (perhaps unexpected) subtleties [29]. Thus, we now consider the 
Lagrangian 

C = V^g[\R - \{dct>f + V,e-'^% (56) 

The resulting equations of motion (the one-field versions of (17) and (18)) 
are solved by the one- field scaling solution (20). It is convenient to define the 
parameter e via 

e^\{l + w). (57) 

In inflation this is the slow-roll parameter; here, since the potential is steep, 
w ^ 1 and e rather represents "fast-roll" . In the scaling solution, it is simply 
given by e = It is also useful to continue the following analysis in terms 
of conformal time r, defined via dt = adr. We will denote derivatives w.r.t 
conformal time by a prime, i.e. ^ =' . In terms of conformal time, the scaling 
solution (20) reads (up to unimportant re-scalings) 

«^(_^)./(i-.), 0,^-^ln(-r), P = E|- (58) 



The single- field scalar perturbation theory was first discussed in [48,49]. In the 
next few paragraphs, we will closely follow the presentation of [18]. Including 
scalar perturbations, the metric can be written as 
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dsVo^ = - (1 + 2Ay)dr2 - 2BYidTdx' 

+ [(1 + 2HLY)5ij + 2HTYij] dx'dx^ 

and the perturbed scalar field as 

= 0o(r)+50(T)y 



(59) 



(60) 



where Y{x), Yi{x), and Yij{x) are scalar harmonics. Here we work in Fourier 

— * 

space, so that we are not writing out the implicit subscripts k of perturbation 
modes. The perturbed Einstein equations SG^^ — ST^^^, then relate the metric 
and the matter perturbations to each other. As is well known, scalar perturba- 
tions in a spatially-fiat FRW universe with scalar field and potential V{(j)) are 
completely characterized by a single gauge- invariant variable. Various choices 
for this variable are typically considered, including the "Newtonian potential" 
and the "curvature perturbation" C,. The gauge-invariant Newtonian po- 
tential $ is most easily understood in "Newtonian gauge" {B — Ht = 0), where 
it is related to the metric perturbations in a simple way: ^ — A — —H^. It 
obeys the equation of motion 



$" + 2 




$' + 2 


+ 271' 




a 0() 







2«f 

00 



$ = 



(61) 



where k = \k\ is the magnitude of the (comoving) Fourier 3- vector. On the 
other hand, the gauge-invariant perturbation variable ( is most easily under- 
stood in "comoving gauge" {Ht — 5T^^ — 0), where it represents the curvature 
perturbation on spatial hypersurfaces, and is related to the spatial metric per- 
turbation in a simple way: ( = —H^. The condition 6T^^ — also implies that 
50 = in this gauge. ( obeys the equation of motion 



C" + 2-C' + k'C = 0, 

z 



(62) 



where we have defined z = a^ct)^/ a' . $ and C are related to each other by 



e 
a' 



— $' + $ 
a' 



ak"^ 



(63) 
(64) 



It is convenient to introduce new variables, u and v [81] (these are of course 
unrelated to the coordinates used in section 3.2), by multiplying $ and C, by 
/c-independent functions of r: 

u = ^$ V = z(. (65) 

00 
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Note that u and v have the same A;-dependence as $ and C, respectively, and 
therefore the same spectral properties, u and v obey the simple equations of 
motion 



^" + (^'-%Tt)^ = (66) 



z" 



v" + {k'' - —)v = Q, (67) 
z 

and are related to each other by 



z' 

kv = 2k(u' + -u) (68) 

z 

When e is time- independent, we can use (20) to find 

(i/zY 



l/z (e-l)V 
z" 2 - e 



1)V2 



(70) 
(71) 



Thus, at early times (when |r| is large) the k'^ term dominates the (l/z)" /{l/z) 
and z" /z terms in the brackets in the equations of motion (66) and (67). Phys- 
ically, this means that the relevant modes are well within the horizon, where 
the solutions are asymptotically oscillatory. On scales much smaller than the 
horizon, the curvature of spacetime is negligible, and thus we impose boundary 
conditions corresponding to the Minkowski vacuum seen by a coming observer 
[14]: 



V -> -^e-'''^ (73) 
V2k 



as r — >• — oo. 



The equations of motion are Bessel equations, and they can be solved exactly 

by 

u(x) = x'/^ [A^'^H^^\x) + A^'^H(^\x)] (74) 
v(x) = x^/^ [b^^^H^^\x) + B^^^Hfix)] (75) 
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where a; = A;|r| is a dimensionless time variable, A^^''^^ and B^^'^^ are constants, 
{x) are Hankel functions, and we have defined 



a = 



e + 1 



e- 1 

e-3 



(76) 
(77) 



In the far past {x — > oo) we use the asymptotic Hankel expression 



exp 

V nx 



, . , sir TT 



(78) 



so that the boundary conditions pick out the solutions 



with the phase factors 



(79) 
(80) 



Vi = exp[i(2a + 3)7r/4] (81) 
P2 =exp[i(2/? + l)7r/4]. (82) 

In order to determine the power spectra of C, and at late times when co- 
moving scales are well outside Hubble radius, we need the asymptotic form of 
the Hankel functions as x — > : 

H^^){x)^--T{s)C-V (83) 

TT Z 

where s > and r(s) is the Euler gamma function. On large scales, the power 
spectrum of the Newtonian potential $ is given by 



^'.W = |^^«-'-" (84) 



The spectral index for $ is thus 

n$-l = l-2Q; = l-|^^|. (85) 

Note that it is invariant under e 1/e, which is a manifestation (at the 
hnear level) of a duality between inflation (e ~ 0) and ekpyrosis (e 1) 
[18]. In fact, this could have been anticipated on the basis of the invariance 
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of both the equation of motion and the boundary condition for u under the 
transformation e 1/e. Note that in both hmits of interest, namely e — 
(inflation) and e — > oo (ekpyrosis), the spectrum of the Newtonian potential 
turns out to be scale-invariant. If we had performed the same calculation, but 
with an equation of state that is slowly varying in time, then we would have 
obtained the spectral index [39] 

n<& - 1 = -4(e + r;), (86) 

where 

e-=(f)^ 9=1-^. (87) 



On large scales, the power spectrum for ( is given by 



nW = ^^ (88) 

which means that the spectral index for ( is 

nc- 1 = 3-2/3 = 3- |^^|. (89) 

In ekpyrotic models, e is large, and thus the spectrum of the curvature per- 
turbation C is very blue, meaning that there is much more power on smaller 
scales, in disagreement with observations. 

This seems very puzzling at first, since one would have thought that it didn't 
matter whether one used $ or ^ in order to perform the calculation. However, 
it is instructive to rewrite the relationship (63) between the two variables in 
the following form: 

1 

c-—,{^y. (90) 

ea^ a' 

Expanding the Hankel function H^^\x) again as x — > 0, but this time up to 
second order, one finds 

Hi'\x) = -;^r(5)(|)-- r(-5)(|)^ + 0(x^-, x'+^), (91) 

which means that $ behaves roughly speaking as 

$ ^ A;-3/2-iA(_^)-i-2A + (92) 

where we have expanded all exponents up to linear order in 1/e. But using 
the same expansion, we have that 



a3 



- (-r)-^-2/^ (93) 
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so that from (90) we can see that the leading, scale-invariant contribution to 
$ drops out of (\ 

Here we see very clearly the importance of adding gravity to our analysis. In 
light of the discussion in section (3.2) this means that the single-field ekpyrotic 
model does not lead to a scale-invariant spectrum of cosmological perturba- 
tions unless there is a mixing of $ and ^ at the bounce (in that case, the 
scale-invariant component of $ would dominate on large scales over the blue 
intrinsic spectrum of Q. Such a mixing would depend on the precise dynamics 
at the bounce, and would therefore be a model-dependent prediction. Exam- 
ples of mixing due to higher-dimensional effects have been given in [78,99], 
and mixing can occur for example due to the breakdown of the 4-dimensional 
effective field theory. But what this means physically is that new degrees of 
freedom become relevant near the brane collision, and therefore it is rather 
natural to try to employ a more complete effective theory, containing more 
fields, in order to perform the above calculation. This will be the topic of the 
following section. 

4.1.2 Two Fields 

In a sense it is rather unnatural to consider only a single scalar field in the 
effective theory, since there are at least two "universal" fields that are always 
present in a higher- dimensional context: the radion field, determining the dis- 
tance between the two end-of-t he- world branes, and the volume modulus of the 
internal manifold. But as soon as there are more than one scalar field present, 
one can have entropy, or isocurvature, perturbations, which are growing mode 
perturbations in a collapsing universe [83] . Entropy perturbations can source 
the curvature perturbation, and hence (provided the entropy perturbations 
acquire a nearly scale- invariant spectrum), almost scale-free curvature pertur- 
bations can be generated just before the bounce. These then turn into growing 
mode perturbations in the ensuing expanding phase. 

As will be shown in section 6, the effective theory arising from heterotic M- 
theory, which describes gravity and the two universal scalars mentioned above, 
is a 4-dimensional theory of gravity minimally coupled to two scalar fields 
(after a field redefinition). It is thus surprisingly simple. On top of that, we 
assume an attractive force between the boundary branes, modeled by two 
negative nearly exponential potentials for the scalar fields: 

V^-Vie-I ^^"^^^ - V2e- / "^^'^^ , (94) 

where Ci = Ci{(f)i), C2 = £2(^2), and Vi and V2 are positive constants. We 
consider potentials in which the q are slowly varying and hence the potentials 
are locally exponential in form. Furthermore, for simplicity, we focus on scaling 
background solutions in which both fields simultaneously diverge to —00. The 
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background solution we arc interested in is the two-field version of the scaling 
solution (20). We will find it useful to define the variable a via [38] 



d = + 0i- (95) 

a has the interpretation of being the path length along the background scalar 
field trajectory. The angle Q of the background trajectory is defined by [38] 
cos(^^) = ^i/o" , sin(^^) = Then the fast-roll parameter e and the equation 

of state parameter w are related to the background evolution via 

= (96) 

During the phase in which the entropic perturbations are generated, ^ = 0, 
which corresponds to a straight background trajectory in scalar field space. In 
this case, we define 

02 = 701, (97) 
and, with this notation, c\ — 7C2 and 

= 2(TT7) ■ ^''^ 

Thus we have an ekpyrotic phase as long as C\C2 is large. 7 is typically of 
C(l); for the heterotic M-theory colliding branes solution with empty branes, 
7 = — ^ during the ekpyrotic phase, as will be derived in section 6. It turns 
out though that the results are not very sensitive to 7. 

The angle of the background trajectory is related to the potential via 

01 +02 

Hence, since we want a straight line trajectory, we have 

T^ = V^(0i)+7'^"(02/7), (100) 

for some function V^(0). 

In a contracting universe, a growing mode is given by the entropy perturbation 
5s, namely the relative fiuctuation in the two fields, defined (at linear order) 
as follows [38] 

5S= (0l502-02 50i)/(^. (101) 

The entropy perturbation is gauge- invariant and it represents the perturbation 
orthogonal to the background scalar field trajectory, as shown in Fig. 6; see 
Ref. [58] for its definition to all orders. 
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The entropy perturbation equation of motion reads (see e.g. [38]) 

5s + ?>H5s + — + + 3^M (^s = , (102) 

Successive derivatives of the potential with respect to the entropy field are 
given by 



HK^. - 02v:<^j (103) 
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(104) 
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+30102 K'/'101'/'2 -02 K'/'Kif-K/-!)- (105) 

Incidentally, note that this implies 9 — —Vg/ a. 

Again, for simplicity wc will focus attention on straight line trajectories in 
scalar field space. Since = 0, the entropy perturbation is not sourccd by the 
Newtonian potential $ and we can solve the equations rather simply. On large 
scales, the linear entropy equation then reduces to 

5s + 2,H5s + Vss5s = 0. (106) 

It is convenient at this point to continue the analysis in terms of conformal 
time T. Introducing the re-scaled entropy field 

5S = a(r) 5s, (107) 

Eq. (106) becomes 

5S" + ^-^ + a^V^,^,^ 5S = 0. (108) 
The crucial term governing the spectrum of the perturbations is then 

2 f ^ T/ 2 



T'^--V^,^,a'y (109) 

When this quantity is approximately 2, we will get nearly scale-invariant per- 
turbations. We proceed [61] by evaluating the quantity in (109) in an ex- 
pansion in inverse powers of e and its derivatives with respect to A^, where 
= hi{a / ttend) 1 where a^nd is the value of a at the end of the ekpyrotic phase. 
Note that A^ decreases as the fields roll downhill and the contracting ekpyrotic 
phase proceeds. 
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We obtain the first term in (109) by differentiating (19), obtaining 

- = 2i/Vfl-^ey (110) 



The second term in (109) is found by differentiating (96) twice with respect to 
time and using the background equations and the definition of N. We obtain 

aV^,^, = -a'H' l2e' - 6e - ^ e,^) + ©(e^^). (Ill) 

Finally, need to express 7i = {o! /a) — aH in terms of the conformal time r. 
From (110) we obtain 

n' = 'H'{l-e), (112) 

which integrates to 

7^-1= r drie-l). (113) 
Jo 

Now, inserting 1 = d(T)/dT under the integral and using integration by parts 
we can re-write this as 

= er (^1 - i - (er)-^ j^^ eVdr) . (114) 

Using the same procedure once more, the integral in this expression can be 
written as 

(er)-! r dr^T = — - (eT)-^ ^(eV) rdr. (115) 
JO e Jo dr 

Now using the fact that e' = Tie^Ar, and that to leading order in 1/e, Ji can be 

replaced by its value in the scaling solution (with constant e), Tir = e~^, we 

can re-write the second term on the right-hand side as 



-7^A(eV)rdr = -(6r)-7^dr-(^) , (116) 
JO dr JO e V e j n 



which shows that this term is of order 1/e^ and can thus be neglected. Alto- 
gether we obtain 



H 



-' = /;dx(.-l)«.r(l-i-if). (117) 



Using (110) and (111) with (117) we can calculate the crucial term entering 
the entropy perturbation equation. 



-=(--K«.«^)=2(l-| + 5S^). (118) 



Following the same steps as in the single field section above, it is straightfor- 
ward to see that the deviation from scale-invariance in the spectral index of 
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the entropy perturbation is now given by 

= (119) 

The first term on the right-hand side is a gravitational contribution, which, 
being positive, tends to make the spectrum blue. The second term is a non- 
gravitational contribution, which tends to make the spectrum red. We will 
return to this expression shortly. Before explaining how these entropic per- 
turbations arc naturally converted to curvature perturbations, we must add 
an important comment: the growth of the entropy perturbations is intimately 
linked to an instabihty of the background trajectory [61,100]. Indeed, the back- 
ground trajectory follows a ridge in the potential [57], and in order to achieve 
a sufficiently long ekpyrotic phase, the trajectory has to remain near the crest 
of the ridge for long enough. This places important constraints on the two-field 
models, which we will discuss in section 5. 

We have shown how an approximately scale-invariant spectrum of entropy 
perturbations may be generated by scalar fields in a contracting universe. 
Wc will now discuss how these perturbations may be converted to curvature 
perturbations if the scalar field undergoes a sudden acceleration, and we will 
estimate the curvature perturbation amplitude. At linear order, a non-zero 
entropy perturbation combined with a bending {9 ^ 0) of the background 
trajectory sources the curvature perturbation on large scales and results in a 
linear, gaussian curvature perturbation (see e.g. [38]) 




where we denote the duration of the conversion by At {i.e. At is the time 
during which 6^0) and the sign corresponds to a contracting universe. Thus 
the strength of conversion is proportional to 1/ and this dependence on the 
equation of state during the time of conversion will have repercussions for the 
magnitude of the second order correction as well, as will be discussed in the 
next section. There are many ways in which such a bending of the background 
trajectory can occur. We will now discuss the various possibilities considered 
in the literature, namely conversion after the ekpyrotic phase during kinetic 
energy domination [61], conversion during the ekpyrotic phase [57,56] or during 
the transition to a ghost condensate phase [22,30], and conversion after the 
big bang by modulated preheating [9]. 
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Fig. 6. The trajectory in field space reflects off a boundary at (^2 = 0. The en- 
tropy perturbation, denoted 6s, is orthogonal to the trajectory. The bending causes 
the conversion of entropy modes into adiabatic modes 6a, which are perturbations 
tangential to the trajectory. 

Conversion during kinetic energy domination 

In the original ekpyrotic and cyclic models, the phase dominated by the steep, 
ekpyrotic potential V{(j)) comes to an end (at t = tend < 0) before the big 
crunch/big bang transition (t = 0), and the universe becomes dominated by 
the kinetic energy of the scalar fields, as we saw in section 3. Consequently, 
the equation of state at t = tend changes from egfc S> 1 to e = 3 (corresponding 
to w — > 1, the equation of state for a kinetic energy dominated universe). The 
conversion from entropic to curvature perturbations then takes place during 
this kinetic energy dominated phase [61]. The conversion occurs naturally in 
the heterotic M-theory embedding of the cyclic model because the negative- 
tension brane bounces off a spacetime singularity [67] - creating a bend in 
the trajectory in field space in the 4-dimensional effective theory - before it 
collides with the positive-tension brane (the big crunch/big bang transition). 
This will be explained in detail in section 6. Here, all we need to know is that 
in the 4-dimensional effective description there are two scalar field moduli, 0i 
and 02, living on the half-plane —00 < (pi < 00, —00 < 02 < 0. In other words, 
there is a boundary to moduli space at 02 = 0. Furthermore, the cosmological 
solution of interest is one in which 02 encounters this boundary 02 = 0, and 
reflects off it at time t = tref', see Figure 6. This bending of the trajectory 
automatically induces the conversion of entropy to curvature perturbations. 
Here we will not restrict the analysis to this particular example, but we will 
consider the general situation in which the scalar field trajectory bends in 
a smooth way during the phase of kinetic energy domination following an 
ekpyrotic phase. 

In general, in the presence of N scalar fields with general Kahler metric gij{4>) 
on scalar field space, the equation for the evolution of the curvature pertur- 
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bation on large scales can be rewritten as [61] 



C^f,«^^. (122) 
where the N — 1 entropy perturbations 

= (123) 

are just the components of 50* orthogonal to the background trajectory, and 
the operator D'^ / Dt^ is the geodesic operator on scalar field space. Here, things 
simplify because the scalar field space is fiat, so the metric is gij = 6ij, and 
D/Dt reduces to an ordinary time derivative. Considering only two scalar 
fields, we have 

= -02 5s/4^^l = +01 8s/4^^l (124) 

As a first approximation, we will assume that the scalar field trajectory simply 

refiects off the boundary at 02 = 0. The scalar field trajectory is 02 = — 70i, 
for t < th, and 02 = 70i, for t > tb, with 0i constant and negative in the 
vicinity of the bounce. The bounce leads to a delta function on the right-hand 
side of (122), 

where from the higher-dimensional point of view tref is the time of the bounce 
of the negative-tension brane. As can be readily seen from (122), if the en- 
tropy perturbations already have acquired a scale-invariant spectrum by the 
time tref, then the bounce leads to their instantaneous conversion into cur- 
vature perturbations with precisely the same long wavelength spectrum. We 
can estimate the amplitude of the resulting curvature perturbation by inte- 
grating equation (123) using (125). Since we have assumed the universe is 
kinetic-dominated at this time, H — l/{3t). As pointed out earlier, since the 
entropy perturbation (95) is canonically normalized, its spectrum is given by 
(55) up to non- scale-invariant corrections. This expression only holds as long 
as the ekpyrotic behavior is still underway: the ekpyrotic phase ends at a time 
tend approximately given by \Vend\ = '^/{cHlnd}- A^^er tend, during the kinetic 
phase, the entropy perturbation obeys 5s -\- t~^5s — 0, which has the solu- 
tion 5s — A + B\n{—t). Matching this solution to the growing mode solution 
t~^ in the ekpyrotic phase, one finds that by tref the entropy grows by an 
additional factor of 1 + \n{tend/tref )- Employing the Friedmann equation to 
relate 02 = 70i to H, putting everything together and restoring both the 
Planck constant and the Planck mass, we find for the variance of the spatial 
curvature perturbation in the scale-invariant case. 
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for the perfectly scale-invariant case. Notice that the result depends only log- 
arithmically on tref' the main dependence is on the minimum value of the 
effective potential and the parameter Ci. Observations on the current Hubble 
horizon indicate A^(A;o) ~ 2.4 x lO"'^ with A;o = 0.002 Mpc"^ [52]. I gnorine; the 
logarithm in (126), this requires 



or approximately the GUT scale. This is of course entirely consistent with the 
heterotic M-thcory setting [8]. Note also that this means that the ekpyrotic 
phase ends at a time of about 10^ Planck times before the big crunch. 

In fact, our approximation of a perfect reflection should really be improved 
upon, especially since eventually we want to consider second-order corrections. 
Specific examples [67] show that in the presence of matter on the branes, one 
expects a smooth repulsive potential and a gradual bending of the trajectory. 
We will discuss this in more detail shortly in section 4.1.3. It does turn out 
however that for gradual reflections, and taking into account the evolution of 
the entropy perturbation during the process of conversion, the efficiency of 
conversion is only reduced by a factor of about 5 or so [65], and hence the 
above rough estimate for the amplitude of the perturbations is rather good. 
Here we simply note that because of the evolution of the entropy perturbation 
during the time of conversion, sharper reflections in fact lead to less efficient 
conversions. Turning this the other way around, one can also say that sharper 
conversions require steeper ekpyrotic potentials (larger values of Ci) to yield 
an amplitude of cosmological fluctuations compatible with observations. 

Conversion during the ekpyrotic phase 

In the "new ekpyrotic" models [22,30,57], typically the conversion of entropy 
to curvature perturbations is assumed to take place because the background 
trajectory switches from the two-fleld unstable scaling solution to a late-time 
single-field attractor solution; in other words, the trajectory starts out close 
to the ridge of the two-field potential and then falls off one of the steep sides, 
either because of the initial conditions or due to an additional feature in the 
potential (for example, one field could see a sharp rise in its potential if it 
enters the ghost condensate phase before the second one) . Adding a feature to 
the potential makes little difference to the results as long as the reflection re- 
mains gradual [65], so in fact we will not consider an additional potential here. 
Compared to conversion during the kinetic phase, the entropy perturbation 
falls off by rather less (see section 4.1.3), but the value of Cc — eek is larger. 
Altogether this implies from (121) that the resulting amplitude of curvature 
perturbations obeys roughly the same constraints as those discussed above, 
see equation (127). 

Conversion via modulated preheating 




(127) 
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Finally, it has been proposed [9] that, instead of converting entropy perturba- 
tions into curvature perturbations before the big crunch/big bang transition, 
the conversion could occur during the phase shortly following the bang through 
modulated reheating. The concept is that massive matter fields are produced 
copiously at the brane collision and dominate the energy density immediately 
after the bang. The massive fields are assumed to couple to ordinary mat- 
ter with a strength proportional to h{()s)j so that their decay into ordinary 
matter occurs at slightly different times depending on the value of 6s. In this 
way, the ordinary matter perturbations inherit the entropic perturbation spec- 
trum. Since the conversion happens while Cc ~ 3, the resulting amplitude of 
the perturbations should be of the same magnitude as for conversion during 
the kinetic phase. Note however that a detailed prediction is made difficult by 
the fact that h is an unknown function of Ss. Indeed, the entropy perturbations 
are converted with an efficiency [9] 



3 \h^s\ 



2 h 



(128) 



In the absence of a concrete model predicting the shape of the function h, it 
is therefore impossible to go beyond the present order-of-magnitude estimate. 

Comparing Predictions for the Spectral Index 

Now that we have seen how curvature perturbations with the right amplitude 
can be generated from a multiple-field ekpyrotic phase, we should add some 
comments about the spectral index of the perturbations, since this is one of 
the quantities that can be measured with great precision. As seen above, the 
curvature perturbations inherit their spectral index from the entropy pertur- 
bations, and therefore it is given by (119). In this section, we analyze this 
relation using several techniques and compare the prediction to those for cur- 
vature perturbations in inflation and for a cyclic model in which the single-fleld 
time-delay fluctuations are converted to curvature perturbations before/at the 
bang [61]. 

As a first approach, let us consider the model-independent estimating proce- 
dure used in Ref. [50]. We begin by re-expressing Eq. (119) in terms oiJ\f, the 
number of e-folds before the end of the ekpyrotic phase (where djV = (e— l)dA^ 
and e » 1): 

2 dine , , 

"'-' = 1-1^- <^^»' 

This expression is identical to the case of the Newtonian potential perturba- 
tions derived in [50], except that the first term has the opposite sign. In this 
expression, e{J\f) measures the equation of state during the ekpyrotic phase, 
which decreases from a value much greater than unity to a value of order unity 
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in the last M e-folds. If we estimate e ^ M"', then the spectral tilt is 

- 1 - (130) 

Here we see that the sign of the tilt is sensitive to a. For nearly exponential 
potentials (a ~ 1), the spectral tilt is ^ 1 + 1/J\f pa 1.02, slightly blue, 
because the first term dominates. However, there are well- motivated examples 
(see below) in which the equation of state does not decrease linearly with M . 
We have introduced a to parameterize these cases. If a > 1.14, the spectral tilt 
is red. For example, = 0.97 for o; 2. These examples represent the range 
that can be achieved for the entropically-induced curvature perturbations in 
the simplest models, roughly 

0.97 <ns< 1.02. (131) 



For comparison, if we use the same estimating procedure for the Newtonian 
potential fluctuations in the cyclic model (assuming they converted to curva- 
ture fluctuations before the bounce), we obtain 0.95 < < 0.97. This range 
agrees with the estimate obtained by an independent analysis based on study- 
ing infiaton potentials directly [19]. Furthermore, as shown in Ref. [50], the 
same range is obtained for time-delay (Newtonian potential) perturbations in 
the cyclic model, due to the duality in the linear perturbation equations dis- 
cussed above [18]. Hence, all estimates are consistent with one another, and we 
can conclude that the range of spectral tilt obtained from entropically-induced 
curvature perturbations is typically bluer by a few percent. 

We note that it is possible that curvature perturbations arc created both by 
the entropic mechanism and by converting Newtonian potential perturbations 
into curvature perturbations through higher-dimensional effects. In this case, 
the cosmologically relevant contribution is the one with the bigger amplitude. 
In particular, the conversion of Newtonian potential perturbations is sensitive 
to the brane collision velocity [99,78], whereas the entropic mechanism is not. 
So, conceivably, either contribution could dominate. 

A second way of analyzing the spectral tilt is to assume a form for the scalar 
field potential. Consider the case where the two fields have steep potentials 
that can be modeled as y(0i) = -Vbe"/'^'^'^ and 02 = 701- Then Eq. (119) 
becomes 

4(1 + 7^) Acs , , 

where we have used the fact that c(0) has the dimensions of inverse mass 
and restored the factors of Planck mass. The presence of Mpi clearly indicates 
that the first term on the right is a gravitational term. It is also the piece 
that makes a blue contribution to the spectral tilt. The second term is a non- 
gravitational term. For a pure exponential potential, which has = 0, the 
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non-gravitational contribution is zero, and the spectrum is slightly blue, as 
our model-independent analysis suggested. For plausible values of c = 20 and 
7 = 1/2, say, the gravitational piece is about one percent and the spectral 
tilt is Us ~ 1.01, also consistent with our earher estimate. However, this case 
with precisely equal to zero is unrealistic. In the cyclic model, for example, 
the steepness of the potential must decrease as the field rolls downhill in order 
that the ekpyrotic phase comes to an end, which corresponds to c,^ > 0. If c(0) 
changes from some initial value c ^ 1 to some value of order unity at the end 
of the ekpyrotic phase after changes by an amount A0, then c^^ ~ c/A0. 
When c is large, the non-gravitational term in Eq. (132) typically dominates 
and the spectral tilt is a few per cent towards the red. For example, suppose 
c oc 0^ and / c(0) d0 ~ 125; then, the spectral tilt is 

- 1 = -0.03^, (133) 

which corresponds to 0.97 < < 1 for positive < /3 < oo, in agreement 
with our earlier estimate. 

Finally, we can estimate the running of the spectral index, given by dn^/dlnA; 
[54]. Again, we can use the model- independent estimating technique used 
above. Then, starting with equation (130) and using J\f — \Q.{aendHend/ O'H) as 
well as A; = aH, we obtain 

dus _ dris _ 2a 

Since we have a > 1, we can see that the running of the spectral index is of 
order 0{M'''^), which is typically ^ 10~^ or less. The current WMAP bounds 
on the running are |dns/dlnA;| < 0.03 [52] and are thus easily satisfied. 



4.. 1.3 Non-Gaussianity 

Cosmological measurements, as performed for example by the WMAP satel- 
lite, are just now becoming sensitive enough to probe the second order correc- 
tions to the background evolution of the universe. If we stop at linear order 
in perturbation theory, then the equations of motion can be derived from a 
quadratic action S. Hence, the fluctuation modes obey gaussian statistics (as 
the probability ~ e^'^). In this sense the second order corrections give us a 
measure of non-gaussianity in the distribution of cosmological perturbations, 
and, in conjunction with precision measurements of the spectral index, mea- 
sures of non-gaussianity provide very sensitive tests of cosmological models. As 
we will see, ekpyrotic models predict a substantial amount of non-gaussianity, 
in contrast to simple single-field inflationary models. This is principally due 
to the fact that the potential is steep, and hence self- interactions of the ekpy- 
rotic field are important, whereas in infiation, the potential is very fiat and 
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self-interactions of the infiaton are negligible. The fact that non-gaussianity 
will be measured with high precision in the near future (with the upcoming 
Planck satellite) and that hints of it have already been seen [110] justifies a 
detailed examination of the ekpyrotic predictions. For conversion during the 
ekpyrotic phase, the order-of-magnitude of non-gaussianity has been estimated 
in [30,23] and more detailed calculations were performed in [55,24]; for con- 
version during the kinetic phase, the non-gaussianity was calculated in [64]. 
The various models have been compared in [65]. 

What we have to do is rather clear: namely, we have to repeat the analysis 
of the previous paragraphs, but this time we need to include all the second 
order terms as well. To this effect, we will decompose the entropy perturbation 
into a linear, gaussian part and a second-order perturbation by writing 6s ~ 
Ss^^^ + Ss^'^\ Its equation of motion, on large scales and up to second order in 
field perturbations is then given by [58] 



The last term in equation (135) is a non-local term proportional to the dif- 
ference in spatial gradients between the linear entropy perturbation and its 
time derivative. This difference evolves as a^'^ [58], so that it remains approx- 
imately constant during the ekpyrotic phase when a is very slowly varying. 
This ends up being exponentially suppressed compared to the entropy pertur- 
bation itself, which grows by a factor of 10^° or more during this same period 
[61]. After the ekpyrotic phase has ended, the non-local term grows because 
a oc (— 1)~^/'^, but this growth is negligible compared to the exponential sup- 
pression during the ekpyrotic phase. Hence, the non-local term can be safely 
neglected. Note that this is a closed equation for the entropy perturbation. 
Here Va denotes a derivative of the potential along the background trajectory. 

It is useful to recast the evolution in terms of the adiabatic and entropic 
variables a and s. Up to unimportant additive constants (that we will fix 
below), they can be defined by 





_ ^HO] 6s('^5s^'^ 
a 2 J 

-] (5.«)^ + -5e(^) = 



(135) 



0101 + 0202 



(136) 



a 
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0102 - 0201 



(137) 
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Then we can expand the potential up to third order as follows [23]: 

= -Voe^^'^il + es' + ^e'l's\ (138) 

where ^3 is of 0(1) for typical potentials (the case of exact exponentials cor- 
responds to = —4-^2/3). The scaling solution (20) can be rewritten as 

a{t) = {-tf' a = -^In (— \/e^^) « = 0- (139) 

The intrinsic non-gaussianity in the entropy perturbation is produced during 
the ekpyrotic phase and can be determined from the equation of motion for 
the entropy field (135), which reduces to 

Ss + 3Hds + V.sSs + WsssiSs'^^^f = 0. (140) 

The last term is 

]^sss = -^V^ (141) 

and thus the solution, at long wavelengths and up to second order in field 
perturbations, is given by [55,64] 



= SSenJ-^ + c{5Send-^f, (142) 

where 

c=^. (143) 

Thus, the intrinsic non-gaussianity present in ekpyrotic models (second term in 
(142) above) is of order 0{.^/e^) and it is due to the steepness of the potentials 
and the resulting self-interactions of the scalar fields. This is in sharp contrast 
with infiation, where the intrinsic non-gaussianity is extremely small due to 
the fiatness of the potential (or, equivalently, e^^/ <^ 1). 

What is measured is not directly the non-gaussianity present in the entropy 
perturbation, but the non-gaussianity it imprints on the curvature pertur- 
bation. Thus, it is important to know the strength with which this intrinsic 

non-gaussianity gets transferred to the curvature perturbation. The time evo- 
lution of the curvature perturbation to second order in field perturbations and 
at long wavelengths is given in FRW time by 



95 s + ^-^[{Vss + 4e2)(5s«)2 - -^5s5sl (144) 
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where we have omitted a non-local term which can be neglected in ekpyrotic 
models for the same reason as the non-local term omitted from (135) above. 
As we saw in equation (121), at linear order, a non-zero entropy perturba- 
tion combined with a bending {9 ^ 0) of the background trajectory sources 
the curvature perturbation on large scales and results in a linear, gaussian 
curvature perturbation C,^ — /a* ■sJj-05s^^\ Thus the strength of conversion is 
proportional to 1/y^ and this dependence on the equation of state signifi- 
cantly influences the magnitude of the second order correction as well, as we 
will see presently. 

As in most inflationary models, the fluctuations are generated by scalar fields 
with canonical kinetic energy density, so they generate non-gaussianity of 
the "local" type, as considered in Refs. [106,53,75,7]. The local wavelength- 
independent non-gaussian contribution to C, can then be characterized in terms 
of the leading linear, gaussian curvature perturbation C,l according to 

C = Cl + \fNLCh (145) 

using the sign convention for wavelength-independent non-gaussianity param- 
eter f^L in [53]. Current observational constraints on f^L have been reported 
by the WMAP team [52] to lie in the range 

-9 < Jnl < 111 (146) 



at the 2(7 level. 



For the two-field models that we are interested in here, Eq. (144) implies that 
the contributions to Jnl can be divided into three parts [64] 



5 r 2H 

f^ntr^ns^c^_ —Qds^''^ (147) 

3^£ JAt (7 

6C,T -J At 0""^ a 

5 



fNL''''''''^^{Ss^'\tend)r, (149) 



6C 



with the total f^L being the sum of all three contributions. J™'^*'^**^ aj^jggg 
from the direct translation of the intrinsic non-linearity present in the entropy 
perturbation into a corresponding non-linearity in the curvature perturbation. 
By contrast, and f^^srated ^^^y |-,g non-zero even if the entropy per- 

turbation were exactly gaussian. Both contributions are due to the non-linear 
relationship between the curvature perturbation and the entropy perturbation, 
as expressed in equation (144) - the difference is that j'^^^^srated ^^^^ generated 
during the ekpyrotic phase and j^^/'^'^**"" during the conversion. 
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The intrinsic non-gaussianity can be estimated by combining equations (121), 
(142), (143) and (147): 

3Cl V Cc 

~ V^ceefe- (150) 

It is given by the geometric mean of the e parameters during the phases 
of generation and conversion of the perturbations. This is perhaps the most 
important result because it fixes a rough magnitude for \fNL\ which can only 
be significantly reduced through accidental cancelations due to the other terms 
or other effects not included in the present analysis. It explains in a nutshell 
why the non-gaussianity in ekpyrotic/cyclic models is necessarily more than 
an order of magnitude greater than in simple inflation models, and it explains 
why the equation of state during conversion can have a significant quantitative 
effect on the prediction. 

To go beyond this qualitative estimate to a more precise one, we need to take 
account of the details of the conversion mechanism. We will again discuss the 
various possibilities considered above in turn. 

Conversion during kinetic energy domination 

We can immediately perform an order-of-magnitude estimate of the non- 
gaussianity generated by this conversion mechanism (described in section 
4.1.2) since during the kinetic phase = 3, and so (150) would lead us to 
expect f^iL to be of order 

/iVL~V^~C»(ci). (151) 



We will treat the bend as if only one field reflects (02), and we will consider the 
case where > 0, see Fig. 6. Other cases can be related to these representative 
examples by changing the coordinate system in field space appropriately. In 
the heterotic M-theory example, the reflection occurs because 02 comes close 
to a boundary of moduh space (02 = 0) and is forced to bounce [67]. For the 
purposes of this study, we will treat the reflection as being due to a potential, 
^^(02); an additional contribution unrelated to the exponential potentials 
that were dominant during the ekpyrotic phase (but which are negligible dur- 
ing the kinetic energy dominated phase). 

It is important to know the evolution of the entropy perturbation during the 
process of conversion. If there is a phase of pure kinetic energy domination 
before the conversion, then the background scalar fleld trajectory is also a 
straight line during this phase, but with the potentials being irrelevant. The 
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equation of motion reduces to 

Ss + -Ss = 0, (152) 

and by matching onto the ekpyrotic solution and its first time derivative at 
t = tend we find 



5s{t)^5s^^\t) + 5s^^\t) 

= SsU^ + In ^) + cSsl,{l + 2 In (153) 

Incidentally, note that 

during the kinetic phase, while during the ekpyrotic phase tSs ~ Ss. This 
observation will simplify our analysis later on. 

During the conversion, even though the kinetic energy of the scalar fields is 
still the dominant contribution to the total energy, the potential V'^(02) that 
causes the bending has a significant influence on the evolution of the entropy 
perturbation. To analyze this, we will approximate (as in [24]) the bending of 
the trajectory to be gradual by taking 9 constant and non-zero for a period of 
time At, starting from t = tref- Note that, assuming the total angle of bending 
is 0{1) radian, we have \6\ ^ l/\tref \ in this case. Then one can relate the 
derivatives of the potential to expressions involving 9, for example 

Assuming a gradual conversion {At ~ tref), we can ignore higher derivatives 
of 9. (In [64] it was shown that sharp transitions lead to unacceptably large 
values of J'nl', hence these cases are of less phenomenological interest.) To 
satisfy the constraint on the amplitude of the curvature perturbation obtained 
by the WMAP observations, we set tref ~ — 10^ Mp^ [61]. At linear order, the 
equation of motion (135) then reads 

5s^'^ + 3H5s^'^ + {9^ + —9)5s^^^ = 0, (156) 

As indicated by equation (154), we can set 5s*' "* = as a first approximation 
and, thus, neglect the damping term in the equation of motion. Also, we will 
simply evaluate the coefficient of the last term midway through the reflection, 
and define 



uj = 9. 



, 1 + . (157) 
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Fig. 7. The evolution of the hnear entropy perturbation during conversion: the 
sohd (blue) line shows the actual evolution calculated numerically, while the dashed 
(purple) line shows the approximate solution (158) with to = 3, 6 = —1/tref and 
t^^f = -400Mp/. For the purposes of illustration, 6s^^\tref) has been normalized 
to 1. 

For a gradual reflection uj ~ (2 — 3)6'. Then the solution for the linear entropy 
perturbation is 

5s(l) = 5s^^\tref) COSw(t - Uef)- (158) 

Instead of continuing to grow logarithmically, the entropy perturbation actu- 
ally falls off during the conversion, see Fig. 7. This has the consequence that 
the conversion from entropy to curvature perturbations is less efficient than 
one might have naively thought. From (121), we can estimate 

= /I^ L ^'^'^ ^ ^i^-'^'^^^'-'f^ ^^^^^ 

where we have used ec ~ 3, which is a good approximation for subdominant 
reflections (and acceptable for the estimating purposes for dominant ones). 
Since uAt ~ 3, the entropy perturbation evolves over nearly a half-cycle and 
consequently sin u; At is a small factor which we will take to be about 1/3 in 
our estimates, which fits well with numerical results. 

We also need to know the evolution of the second order entropy perturbation 
during the time of conversion. Using the same approximations as above (im- 
plying that we can neglect Vgss compared to OVss/o' and 6^/& by the time the 
conversion is underway), the equation of motion (135) simplifies to 

Ss^'^+M^^-^u;\6s^'^Y = 0, (160) 
a 

where Ss^^^ is given in (158). Putting 6 ^ a at the start of the reflection and 

• (2) 

keeping in mind that we impose the boundary condition 6s ~ 0, the solution 
for the second order entropy perturbation is given by 
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+ ^('^^^'^(W))'( - 4cos[a;(t - Uef)] 
+4cos^[a;(t — tre/)] + 9sin^[ci;(t — tre/)] 

+ sin[a;(i — Uef)] sin[3a;(i — tref)]j ■ (161) 

At large eek (or, equivalently, large |c|), the second order entropy perturbation 
falls off in the same way as the linear perturbation, but at small |c| there are 
significant corrections to this behavior. We will also need the integral 



/■ J„(2) _ (1 + 2 sill(^'At) ^ 
J At {l+ln{tend/tref))'^ ^ 

At sm{ujAt) sm{2ujAt) 

where C = Ss^^\tref)- In all cases of interest, the last two terms are neghgible. 

We arc finally in a position to evaluate the various contributions to the non- 
linearity parameter fNL. The intrinsic contribution, defined in (147), becomes 



fintrinsic^^^^^^^j^^ (163) 

with 

5cj(l + 21n(tend/W)) 



8^6^(1 + In {tend/tref)y Sm{ojAt) 

B^^J^, • (164) 

Eq. (163) is one of our key results because it shows that the essential con- 
tribution scales in a simple way with e^k and has a value that exceeds the 
total /atl for simple inflationary models by more than an order of magnitude. 
As suggested by Eq. (150), the first term in fj^*^^^^^<^ can be re-expressed as 
C {-s/^c^ek) with €c — S. This contribution comes directly from the non-zero 
generated during the ekpyrotic phase, which is due to the third derivative 
w.r.t. s of the ekpyrotic potential. It is interesting to note that, in the case 
where Vgss = or ^3 = 0, jg nevertheless non-negligible because of 

the positive offset B generated by the linear entropy perturbation Ss^^^ during 
conversion, the piece proportional to {Ss^^^Y in (161). 

Using the deflnition (148) together with (155) it is straightforward to see that 
for conversion during the kinetic phase, is always negative, and it 

can be estimated as 
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5 

15cj2 \t^^j. + At/2| 



(165) 



where we have used ■jtrefd ~ 1 and 

/ tsm^u{t-tref)^Utref + ^)At, (166) 

which is a vaUd approximation since 6s^^^ evolves over approximately a half- 
cycle. We have neglected the contribution due to the term proportional to 
dsds, since it is suppressed on account of (154). The integrated contribution 
to fwL generated during the ekpyrotic phase gives an additional contribution 
of 



(167) 



8^2(1 + In {tend/tref)y Sm^CuAt) 



Note that neither the reflected nor the integrated contributions depend on eg^; 
they both simply shift the final result by a number depending on the sharpness 
of the transition and the duration of the purely kinetic phase respectively. 

The total Jnl is the sum of all the above contributions. Since we are only 
considering gradual conversions, we expect the ratio \tref + At/2 1 /At to lie 
between 1 and 2. Also, the duration of the pure kinetic phase between the end 
of the ekpyrotic phase and the conversion is necessarily rather short, so that 
we expect ln(tend/W) ~ ^(1)- 

Putting everything together, we get the following estimates: 



fntrinsic^-^^Q-^^Q (168) 
j^re f lection ^ (169) 
jintegrated^^^ (170) 



The analytic estimate for f^j^l suggests a range that could extend to —100, 
but the value above is more consistent with the exact numerical calculations 
for a wide range of parameters. Altogether, we end up with the following fitting 
formula 
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Fig. 8. A comparison of the results from numerical calculations with the fitting 
formula given in Eq. (172) and indicated by the thick black line. Here we have 
fixed the value of e^k = 36. The plot confirms that /tvl then grows linearly with 
^3, in good agreement with (172). The sample models (dashed and dotted lines) 
are representative of the range of models and parameters shown in [64]. We have 
similarly checked the dependence on y/e when the value of K3 is kept fixed. 

^ 10c + 5 (171) 

3 , , , 

~-K3V^ + 5. (172) 

Fig. 8 shows that this formula agrees well with the results from exact numerical 
calculations. In particular, we have been able to estimate the slope, and hence 
the dependence on KsY^Cefc, rather accurately. 



Conversion during the ekpyrotic phase 

If the conversion happens during the ekpyrotic phase, i.e. while the ekpyrotic 
potentials are relevant, ec = eek, and from (150) we expect the non-gaussianity 
to be of order 

fNL ~ ~ 0{cl), (173) 

i. e. we expect the non-gaussianity to be larger by a factor of ci compared to 
the case of conversion during the kinetic phase. 

As shown by Koyama et al. [55], the formalism is well suited to treating 
this case, and the derivation surprisingly quick. According to the 5M formalism 
[91,87], the comoving curvature perturbation is identified with a perturbation 
in the number of e-folds of ekpyrosis (the perturbations are identified on a 
uniform density hypersurface at t = tf, starting the evolution from a flat 
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initial hypersurface at t — ti) 



C^SAf^ + \Ka.aS^O^g)\ (174) 

where we have expanded SAf up to second order in terms of a gaussian variable 
ag. We will specify the meaning of ctg shortly. Then, comparing with the 
definition (145), we immediately see that the non-gaussianity parameter is 
given by 

We can evaluate N by using the one and two-field versions of the scaling 
solution (20), with the result that [55] 

jV = / Hdt + / ^ Hdt = -4 In \Href \ + constant, (176) 

where the index j corresponds to the field that becomes frozen in the late- 
time single-field solution. Here we have assumed that at time t^e/ the trajec- 
tory instantaneously switches from the unstable two-field scaling solution to a 
single-field attractor scaling solution [57]. Now, using that 5s oc t~^ oc if, we 
can rewrite equation (142) as 

5s = agH + c{agHf, (177) 

whence ag is seen to be a parameter distinguishing different trajectories. Turn- 
ing this relation around, we have that 

5s 

ag^-{l-c5s). (178) 

Hence at the time of transition (reficction) tref, when 5s = 5sref is fixed, we 
have that ag oc H~^Y Thus, using (176), we can rewrite (174) as 

Cjag CjUg 

It immediately follows that the non-gaussianity parameter /jvl is given by [55] 

where, again, the index j corresponds to the field that becomes frozen in the 
late-time single-field solution. This expression has been checked numerically, 
both in [55] and in [65], by different methods, and the numerics have been 
found to agree well with the analytic formula above. Note the particular feature 
that the sign of /atl is always positive for these cases. 
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A qualitative understanding of this result can be achieved in the present con- 
text as follows: from (147) we can see that, since H/a is negative and approxi- 
mately constant, the sign of y™^''*"'**^ jg giyg^ by the sign of / 96s^'^^ . In practice, 
numerical simulations indicate that, for conversion during the ekpyrotic phase, 
we have 

1^1 ~ (181) 

' ' 10|t| ^ ' 

during most of the conversion, but with 6 growing faster towards the end. We 
will use this numerical input to guide our analysis. The sign of y^^''*"'**c jg 
essentially determined by the sign of 9Ss^'^^ towards the end of the period of 
conversion. Naively it is difficult to perform a purely analytic estimate of /tvl 
in the current scheme, since all the terms in the equation of motion (135) go 
as Cit~'^; however, equation (181) tells us that 

= 2/^2 > (182) 

and so there are surprisingly few terms in the equation of motion for ^s*-^-* that 
are actually important during (most of) the time of conversion. In fact, to a 
first approximation, we are simply left with 

Ssl^) = (Avsss + ^V,sW^y. (183) 
z a 

Vsss decreases in importance as the single-field scaling solution is reached. 
Thus, even though the Vsss term determines the initial evolution of Ss^^\ 
eventually the Vss term dominates. Then, since V^s < it is easy to see that 
the sign of Ss^'^^ is always driven to be opposite to that of 9, and consequently 
pntrtnstc jg negative in all cases. 

Since jg proportional to Cl^j it provides a contribution of the same 

order 0{€ek)- Using (182), it is straightforward to see that the part of y^^jf'^"^**"" 
that is proportional to (ds)^ is always positive, while the part proportional 
to dsds is always negative. This implies that there will be a competition 
between the two contributions. Numerical integration then shows that the 
part proportional to (ds)^ approximately cancels out while the dsds 

part by itself is very close in numerical value to the final answer. In all cases 

^integrated (.Qj^^plg^gly negligible. 

Consequently, the total /jvl turns out to be negative and moreover in good 
agreement with the 6N result (180) [55], as shown in Table 1. Clearly, the 
5N formalism provides a fast and elegant derivation of non-gaussianity for 
conversion during the ekpyrotic phase that agrees well with direct integration 
of the equations of motion. (Note, however, no analogous 5N approach applies 
to conversion in the kinetic energy dominated phase.) 

The conclusion that /jvl is substantially less than zero differs from Ref. [24]; 
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Table 1 



Ekpyrotic conversion: the values of /nl estimated by the dA'^ formalism compared 
to the numerical results obtained by directly integrating the equations of motion. 

their approximation method was incomplete though, since the term propor- 
tional to dsds in (144) was not included. As shown above, this term typically 
contributes significantly to the final result. Also, the evolution of the entropy 
perturbation during conversion was neglected, which meant, for example, that 
values of y^l""^*'^ of either sign were obtained. Our results here show that 
their approximations were too crude, and that Jnl is generally negative with 
Inl ^ —20. This is inconsistent with current limits by roughly 3 a [52]. 

Hence, wc conclude that models with conversion during the ekpyrotic phase 
(ec = Cefc ^ 1) are difficult to accommodate with current observations in 
contrast to models in which conversion occurs during the kinetic energy dom- 
inated phase (ec = 3). 

Conversion after the crunch/bang transition 

If the conversion of entropy to curvature perturbations occurs via the scenario 
of modulated preheating after the big bang [9], then the conversion happens 
while ec ~ 3, as described in section 4.1.2. Thus we can expect an intrinsic 
contribution to f^L oi 

fNL^V^^Oic), (184) 

i-S- fNL is of the same order as in the case of conversion during the kinetic 
phase preceding the big crunch. Moreover, as shown in [9], no significant ad- 
ditional contributions to /jv^ are expected. Therefore, models of this type are 
subject to roughly the same observational constraints as the models where the 
conversion happens during the kinetic phase before the bang and where the 
intrinsic contribution to the non-gaussianity is dominant. Note however that 
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the entropy perturbations are converted with an unknown efficiency [9] 



e = - 



3 \h^s 



(185) 



2 h 



Since the non-hnearity in the entropy field is of magnitude c, this implies a 
non-linearity in the curvature perturbation given by 



where the ± sign reflects our ignorance of the sign of h g and of the direction 
of bending of the scalar field trajectory. Thus, in the absence of a more de- 
tailed model, we cannot go beyond this rough order-of-magnitude estimate at 
present. 



The analysis of non-gaussianity can be summarized in a few rules of thumb: 

• The intrinsic contribution to I/atlI is proportional to the geometric mean of 
Eek and Cc, which is at least two orders of magnitude greater than in simple 
inflationary models (where _ 0(o.l)). When all contributions are 
considered, the total f^L is generically more than an order of magnitude 
greater than in simple inflationary models, (where /^^"' = 0{1)). 

• For a flxed value of ^3, the value of is correlated with the spectral 
tilt: smaller imphes smaller egfe, which tends to make the spectral 
tilt bluer. Current limits on /^vl fit well with limits on the spectral tilt for 
the simplest models with conversion during the kinetic energy driven phase 
before the bang or during reheating after the big bang. 

• Models in which the conversion occurs during a phase with larger Cc produce 
a larger intrinsic IJnlI and are more difficult to fit with the current observed 
limits on /jvl and spectral tilt. In particular, models with conversion in the 
ekpyrotic phase favor I/tvlI large and negative. 

• Cases in which the intrinsic contribution to /^vl is much smaller than 
the reflection plus integrated contributions produces very large values of 
IInlI 3> 100 that are inconsistent with current observational bounds. This 
includes all cases where the conversion is sharp. 

The analysis suggests a useful characteristic plot for differentiating cosmo- 
logical models: /tv^ versus tilt. Figure 9 illustrates the prediction for a cyclic 
model in which the conversion from entropic to curvature perturbations occurs 
in the kinetic energy dominated epoch following the ekpyrotic phase. Here we 
keep the skewness K3 of the potential flxed, as we vary its steepness e. The 
prediction is a swath whose width is largely due to the uncertainty in the 
trajectory, parameterized by 7. Although the swath includes Jnl near zero, 
positive /iVL between 10 and 100 is preferred for tilts in the range suggested 




(186) 
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Fig. 9. A plot for characterizing the correlation between {/nlI and scalar spectral 
tilt, Us — 1, here illustrated for the case of the cyclic model in which the conversion 
from entropic to curvature perturbations occurs during the kinetic energy dominated 
phase just before the big crunch/big bang transition. Different curves correspond to 
different fixed amounts of skewness K3 in the potential (the central curve corresponds 
to K3 = 4y^2/3), while we vary the steepness of the potential Cek- The curves show 
the general trend that {/nlI increases as the spectrum becomes redder. Simple 
inflationary models correspond to the narrow horizontal hashed (red) strip with 
I/tvlI ^ 1- The shaded rectangle represents the current observational constraints on 
fNL and tilt (95% confidence) from WMAP5 [52]. 



by WMAP5. The prediction for simple inflationary models is confined to the 
narrow band Af^L ^ 1 around zero. 

The results are surprisingly predictive. If observations of Jnl lie in the range 
predicted by the intrinsic contribution of either inflationary or ekpyrotic/ cyclic 
models, it is reasonable to apply Occam's razor and Bayesian analysis to favor 
one cosmological model over the other. Combining with measurements of the 
spectral tilt significantly sharpens the test. The current observational bounds 
obtained by the WMAP satellite are still inconclusive, but it is clear from 
the estimates presented here that non-gaussianity should be detected by the 
Planck satellite if the ekpyrotic/cyclic model is correct. At the same time, this 
provides a strong incentive to further refine other methods of measuring non- 
gaussianity, such as looking for evidence in measurements of the large scale 
structure of the universe [89]. As an example, a positive f^i would mean that 
the largest structures formed slightly earlier, and thus would have drawn more 
mass out of the voids, and hence one might for example expect slightly emptier 
voids, as compared to a perfectly gaussian distribution of primordial density 
fluctuations. 
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4-2 Tensor perturbations 



The ekpyrotic phase also generates gravitational waves, but with a spectrum 
radically different from the inflationary prediction. Thus, gravitational waves 
provide another useful test for differentiating between models of the early 
universe. Unfortunately, it is unclear at present how soon measurement devices 
will be sensitive enough to probe the gravitational wave background. 

Tensor perturbations are simpler to analyze than scalar perturbations since 
there are no issues of gauge variance. In the context of ekpyrotic models, the 
analysis was first performed in [48,17]. Here, we will follow the treatment of 
[18]. The perturbed metric is given by 

dsVa^ = -dr^ + [5ij + 2hTYlp]dx'dx^ (187) 

(2) 

where Y^j is a tensor harmonic. The tensor perturbation Ht is gauge-invariant 
and obeys the equation of motion 

h'i. + 2{a'/a)h'T + k'^hT^0. (188) 

It is useful to define the re-scaled variable /r = ahr which obeys 

fff+{k''-a"/a)fT^O. (189) 

Again, the standard vacuum choice is the Minkowski vacuum of a comoving 
observer in the far past, corresponding to the boundary condition 




(190) 



We can now solve for /t, just as in the scalar case. But it is simpler to notice 
that (58) implies z{t) oc a(T), and hence z"/z = a" /a when e is constant. 
Thus, since v and fx obey identical equations of motion (compare (67) with 
(189)) and boundary conditions (compare (73) with (190)), we find 

h^v^V.sj^Hfix). (191) 



The tensor spectral index is defined in the long- wavelength hmit by /c^|/tP oc 
k^"^ . Using the small x expansion of the Hankel function and (191) we find 



riT = 3 - 2/? = 3 - 



e- 1 



(192) 



Thus, in ekpyrotic models, the tensor spectrum is very blue, rir ~ 2. There 
is a cut-off of the spectrum at high frequencies [17], due to a combination of 
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two facts: first, gravitational waves arc not confined to the brane - they see 
the whole spacetimc, as they arc excitations of the metric. And secondly, af- 
ter the ekpyrotic phase has ended, the kinetic dominated phase starts on the 
brane worldvolume. However, as described in section 3.2, the corresponding 
higher-dimensional spacetime near the brane collision is compactified Milne 
spacetime, which is locally Minkowski space. Thus, gravitational waves with 
wavelengths smaller than about the GUT scale, which leave the horizon after 
the end of the ekpyrotic phase, behave as in Minkowski space and are not am- 
phfied (see also section 6 for more details on the higher-dimensional spacetime 
near the collision). The high-frequency modes just below the cut-off lead to the 
strongest constraint on the gravitational wave spectrum, as they are suscepti- 
ble of interfering with big bang nucleosynthesis. The corresponding constraint 
requires a minimum temperature at reheating, and will be presented in sec- 
tion 5 along with other constraints that ekpyrotic and cyclic models have to 
satisfy. 



The ekpyrotic spectrum is in marked contrast with the model of inflation, 
which predicts a nearly scale-invariant spectrum of gravitational waves, essen- 
tially because there, the gravitational waves are produced in the same way 
as the scalar perturbations. Because of the very blue spectrum, gravity waves 
are just about unobservable in ekpyrotic models [48]: if we assume the same 
Hubble parameter H^, at reheating in a model of inflation, and at the brane 
collision in an ekpyrotic model, then modes of wavelength H^^ will have ap- 
proximately the same amplitude, namely Hb/Mpi. But now consider a mode of 
wavelength comparable to the present-day horizon, ~ Hq^. Then, for inflation, 
we would expect an amplitTidc of similar magnitude due to the approximate 
scale-invariance of the spectrum. But for ekpyrosis, the current horizon is of 
the order of 60 e- folds larger than H^^, and hence, due to the blue spectrum, 
the amplitude of a horizon-sized gravity wave today would be exponentially 
small and thus unobservable. In fact, on large scales, one would expect the 
leading contribution to the gravitational wave background to be induced by 
the backreaction on the geometry of the first-order scalar perturbations [11]. 



Hence gravitational waves provide one of the crucial observational tests that 
are capable of unambiguously ruling out ekpyrotic models. If gravitational 
waves are detected, and it can be established that their spectrum is close to 
scale- invariant, this will constitute very strong evidence in favor of inflation. If, 
on the other hand, primordial gravitational waves remain unobserved, ekpy- 
rotic models will deserve our attention and we should look for other signiflcant 
tests (such as non-gaussianity) in order to improve our understanding of the 
early universe. 
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5 The End is the Beginning: a CycUc Model of the Universe 

We have seen how the ekpyrotic phase can resolve the standard puzzles of big 
bang cosmology, and how it can generate scalar and tensor fluctuations with 
interesting properties. However, this leaves many aspects of our universe un- 
addressed, and foremost the rather mysterious dark energy. The cyclic model 
[92-96,51,35] is an ambitious attempt at providing a complete history of the 
universe, within the framework of a braneworld view of the universe, while 
incorporating both the ekpyrotic mechanism and dark energy in an essential 
way. We will first outline the main ideas, then discuss the constraints that 
the cyclic model has to satisfy and finally describe the global view and new 
possibilities that the model offers. The cyclic model relies on having a scalar 
field potential of the general shape depicted in figure 10. The plausibihty of 
such a potential arising within the context of string theory will be discussed 
in section 6. The idea is that the ekpyrotic phase, which precedes the big 
crunch/big bang transition, is itself preceded by a phase where the potential 
is positive and fiat. During this phase, the radion field rolls down the potential 
very slowly, and this in fact provides a period of dark energy domination (in 
the form of quintessence [84]). But this could be our current universe! The 
dark energy that we see today is reinterpreted as a small attractive force be- 
tween our brane and the second, parallel, end-of-the-world brane. At some 
point in the future, this force becomes stronger; the potential energy becomes 
negative, and the universe reverts from expansion to contraction. We then 
enter an ekpyrotic phase, which locally fiattens and homogenizes the branes. 
After a brief phase of kinetic energy domination, the branes collide, matter 
and radiation are produced and brane dwellers experience a big bang. After 
the bang, the scalar field acquires a small boost and quickly rolls back across 
the potential well and onto the positive energy plateau, where the scalar gets 
quasi-stabilized. Now the universe undergoes the usual periods of radiation 
and matter domination, while the scalar field starts rolling slowly back. Even- 
tually, the resulting dark energy comes to dominate the energy density of 
the universe, and the whole cycle starts again. In this way, physical processes 
in the universe today provide the initial conditions for the next cycle. What 
should be noted is the amazing economy of ingredients: only branes and their 
motion are used to address the general evolution of the universe over its entire 
history! Of course, aesthetic appeal is nice in a physical model, but we must 
ask the question: can this work? In order to answer this, we must have a closer 
look at the details. 

First, there is the question of entropy [93]. The second law of thermodynamics 
implies that the entropy of each cycle should be larger than that of the previous 
cycle; but this seems to preclude the possibility of a truly cyclic universe 
(and historically, Tolman used this reasoning to show the impossibility of 
cychc models in closed universes with zero cosmological constant [101]). The 
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Fig. 10. The potential for the cyclic universe integrates the ekpyrotic part and a 
quintessence epoch, but is irrelevant at the brane collision. A possible form for the 
potential is V{(j)) = Vb(e'"^ - e-''^)F{(j)), with 6 <C 1, O 1 and F(0) tends to unity 
for (p > (pend and to zero for (j) < (kend- Reproduced with permission from [35]. 

resolution is that, as we will see shortly, there is a huge net expansion every 
cycle, and hence, even though the total entropy always increases, the entropy 
gets diluted during the dark energy phase, and the entropy density can return 
to approximately the same value cycle after cycle. Note that the ekpyrotic 
phase does not undo the dilution, since the scale factor shrinks very modestly 
during ekpyrosis. 

It is useful to keep track of various scales during the evolution of a given 
cycle. We are now at the start of a dark energy phase, during which the 
scale factor a grows by an amount e^°^ ^ where, by definition, Mbe is the 
number of e-folds of dark energy domination. During this phase, the Hubble 
parameter remains approximately constant. Then, during the ekpyrotic phase, 
the scale factor shrinks very modestly, by iKnd/Hegl"^^^ ~ ^(1); where V^eg 
is the value of the potential at the beginning of the ekpyrotic phase; this value 
coincides with the value of the present day cosmological constant Vfcej, ~ Vq. 
During the same time, the Hubble parameter grows immensely, by an amount 
I Vend/Hegl^'^^ = e-^<='=, wheie AAefc deuotes the number of e-folds of ekpyrosis. As 
the branes collide, radiation and matter is produced and the branes "reheat" 
with a finite temperature, where we define to be the temperature at the 
time that the radiation density equals the scalar kinetic energy density. There 
is a contracting kinetic energy phase just before the brane collision, followed 
by an expanding kinetic phase after the bang. These two phases are almost 
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exactly symmetrical, except for a slight boost that the scalar field obtains 
at the boimce. This boost must occur in order to achieve a cyclic model, 
because, due to extra Hubble damping from the radiation, the scalar needs a 
slight excess of kinetic energy to cross the potential well onto the plateau, see 
figure 10. There are two known ways by which the scalar field can obtain this 
slight excess energy [93,104]: either due to the coupling of the matter fields 
to the radion, or if there is slightly more matter produced on the negative- 
tension brane than on the positive-tension one. The latter result makes sense 
heuristically, as more negative than positive inertia is produced in that case. 
The first expanding kinetic phase is followed by a brief, and unimportant [35] , 
w ^ 1 phase, which in turn is followed by another expanding kinetic phase. 
We can treat all of these phases as a single kinetic phase, as in [35]. Then, 
since p oc oc during kinetic energy domination, the net expansion of 
the scale factor is given by (| KndT^'^/^r)^^^ = ^ikeI'^^ where we define the 
parameter 

IKE ^ In (^^) ■ (193) 

Meanwhile, the Hubble parameter has shrunk by (| Kndl^/^/Tr)"^ = e~^^^^. 
Finally, during radiation and matter domination, the universe grows by a 
factor of Tf./To = e^™'*, where Tq is the temperature of the CMB today; at the 
same time, the Hubble parameter shrinks by (Tj./To)~^ = e"^^'''"'. Therefore, 
to recapitulate, the scale factor a grows by a total of 

A/'DB+^ + A^.ad (194) 

e-folds over the course of a single cycle. However, because 



C (-Kn.)V^T2 ' 

due to Vfteg ~ Tq , the Hubble parameter actually returns to its original value 
after one cycle. This result can also be re-written as the statement that [35] 

Kk^'^inKE + Nrad). (196) 

There is no constraint on the number of e-folds of dark energy domination 
Nde- 

In order for the cyclic model to work, the temperature has to satisfy a 
number of constraints. The most obvious one is that the reheat temperature 
has to be high enough for nucleosynthesis to take place. Also, it should be 
below the GUT scale, so that no topological defects are formed. Hence the 
first constraint is 

IMeV <Tr< lO^^GeV. (197) 
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Next, in order for the successful predictions of nucleosynthesis not to be ad- 
versely affected by the overproduction of gravitational waves [5,90], the tem- 
perature has to satisfy [17] 

Tr ^ ^\Vendl (198) 

where p is defined in (20) and is typically of (9(10^^) for realistic models. Also, 
the radiation component should not dominate the cosmological evolution until 
the scalar field has crossed the potential well and gotten onto the plateau of 
the potential. Now, combining the evolution of the scalar field, equation (32), 

with the form of the ekpyrotic potential V ~ —e"'^^, we have that V oc f^f^'^ . 
Hence, the time tf^oss it takes for the scalar field to cross the well [i.e. to go 

from — Vaid to Vo) is t cross = ^endlKnd/Vbl'^^^^"^. Wc uccd the time of radiation 
domination tr to be larger than tcross-, or, using that tr ~ H^^ ^ T^"^^ we 
get the constraint (this bound was derived in [51], although there is a typo 
regarding the sign of the inequality in the paper) 

Tr<{-V,ndf"\^\^'. (199) 
Vend 

It should be noted that a broad range of parameters satisfies all the constraints 
for cycling [51]. Typical values for the various parameters are for example 

\V^^^\ = (10^^G'ey)^ Tr = lO^^Gey, p = 10"^ (200) 

Then, with Tq = lO^^eV, we have Jke ^ 7 and N^ad ~ 55. Hence this leads 
to more than 120 e-folds of ekpyrosis, largely sufficient to solve the standard 
cosmological puzzles. And assuming zero e-folds of dark energy domination, 
according to (194), the universe grows by a factor of about ~ 10^^ per 
cycle. Note that this means that our currently observable universe measured 
only about a cubic meter one cycle ago! 

The enormous time scales involved in the cyclic model (a single cycle can al- 
ready easily last a trillion years) mean that new approaches are possible for 
addressing some long-standing problems, such as for example the cosmological 
constant problem. In trying to explain the smallness of the cosmological con- 
stant by a dynamical relaxation mechanism, the problem has typically been 
that there simply isn't enough time between the big bang and now in order 
to achieve the required level of tuning. This has been an issue with the model 
of Abbott [2] for example. However, in the cyclic model, the same model can 
work, as shown in [97]. The idea is that the plateau in the cyclic potential 
is really part of a tilted washboard-shaped potential for an axion-type scalar. 
That scalar can then tunnel down the potential, spending more and more 
time at lower and lower values of the cosmological constant. It spends an ex- 
ponentially longer time at the lowest positive local minimum of the potential 
(before tunneling to negative values, where the universe then collapses), and 
hence that value of the cosmological constant is the most probable one. This 
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is just one example, but the general idea is that there could be evolutionary 
time scales vastly larger than the ones that are usually considered, and which 
can lead to qualitatively new explanations. 

Let us turn to the question of "initial" conditions. The single-field cychc model 
is an attractor, while at the same time managing to generate scalar perturba- 
tions that are nearly scale-invariant [93,42,39]. This is really quite remarkable, 
and it means that the issue of initial conditions is essentially circumvented. 
Evidently, avoiding the issue of initial conditions is one of the main attrac- 
tions of cyclic models. However, for two-field models [61], the situation is more 
complicated. As discussed in section 4.1.2, in order to achieve a nearly scale- 
invariant spectrum of scalar perturbations, the background trajectory must 
fall off along a ridge of the potential, i.e. in the presence of an unstable trans- 
verse direction. This instability is essential to obtain perturbations that are 
compatible with observations [100], and hence there is an issue of initial con- 
ditions. This is an open problem at present. The question of how the field 
returns to a starting position very close to the ridge after each cycle remains 
to be addressed. Of course, to some extent this problem is linked to the gen- 
eral problem of moduli stabilization (or quasi-stabilization in the cosmological 
context). In heterotic string theory, this problem is less well understood than 
for example in the type IIB theory at present. More work is certainly needed 
in this area. In the more phenomenological context of the new ekpyrotic mod- 
els, Buchbinder et al. [23] have proposed to include a stabilizing positive mass 
potential at the top of the ekpyrotic ridge, while coupling the entropy field to 
light fermions. Then, during this pre-ekpyrotic phase, the entropy field oscil- 
lates around the minimum of the potential. The field gets localized near the 
top of the ridge, because the oscillations are damped as energy is dissipated in 
the form of thermal radiation due to the coupling of the entropy field to light 
fermions. Something similar could also work in the cyclic model, and in fact, 
as the entropy field rolls slightly up the stabilizing potential, this increased 
potential energy could be identified with today's cosmological constant ^ . Of 
course, one would like to see such a potential arise from a fundamental theory, 
and avoid putting it in by hand. 

Finally, we will address the global structure of the cychc universe and some 
more speculative questions. Since the universe undergoes an exponentially 
large net expansion every cycle, over large patches the universe will approxi- 
mately have the structure of de Sitter space, but with singular hypersurfaces 
corresponding to the big crunch/big bang transitions [96]. Hence, starting 
from a given patch of spacetime, the Penrose diagram will look approximately 
like the one shown in Fig. 11. The cycling can however not be synchronous 
over super-horizon scales, as the ekpyrotic phase only homogenizes the branes 
over large, but finite-sized patches. Occasionally, there will be large quantum 



^ This idea was suggested to the author by Paul Steinhardt. 
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Fig. 11. On average, the Hubble parameter and the energy density are positive in 
the cychc universe, and hence it can be represented approximately by the conformal 
diagram of de Sitter space, but with singular hypersurfaces corresponding to the 
big crunch/big bang transitions (here these are drawn only in the region covered 
by the usual flat coordinates on de Sitter space). In this diagram, time is vertical, 
and space horizontal, with each point representing a 2-sphere. A representative light 
signal is represented by the dashed line - it travels through many crunches, which 
makes it unlikely that one could see back through many cycles, even in principle. 

fluctuations as the brane collision is approached, leading to the formation of a 
large black hole at the collision. In this region, the cycling will stop [35]; but 
since there is such a huge expansion of the universe during each cycle, these 
regions will be exponentially rare. In the same vein, one might ask whether an 
advanced civilization, with sufficient control over the spacetime metric, could 
stop the cycling in order to prolong its lifetime. Presumably so, but due to 
causality, only the region in its vicinity would be affected, and again, due to 
the vast net expansion, there would be plenty of room left for the cycling to 
continue unaffected. Thus, during each cycle, vast new regions of space are cre- 
ated, but within the same universe; in contrast to eternal inflation [70], there is 
no proliferation of universes. Note also that it seems likely that only informa- 
tion from the previous cycle is remembered, and not from the ones before [35]. 
This is essentially due to the effectiveness of the ekpyrotic phase in washing 
out all classical inhomogeneities, while new inhomogeneities are produced by 
random quantum processes. Nevertheless, could an advanced civilization have 
left us a message from the previous cycle? This seems extremely unlikely, as 
it would need to have had access to precisely the one cubic meter of space 
that has now become our currently observable universe. Could there have 
been a beginning to the cycling, for example if the two branes tunneled into 
existence? This is certainly possible, but if one really cannot see back more 
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than one cycle prior to the current one, then this question might be rather 
meaningless. Finally, we note that, according to the present understanding, 
there seems to be no limit to the number of cycles allowed ^ . Hence the cycles 
can continue indefinitely, creating more and more space, while injecting new 
matter and life into the universe periodically. 



6 The Link to Fundamental Theory: The Embedding in Heterotic 
M-theory 



It is the link between 11-dimensional supergravity and heterotic string theory 
which led to the colliding branes picture in the first place, and even though 
many of the relevant calculations are done in a 4-dimensional effective field 
theory, the higher- dimensional picture is of crucial importance for a deeper 
understanding of ekpyrotic/cyclic cosmology. Moreover, it is helpful in trying 
to visualize many of the processes described in the rest of this review, as it 
provides a geometrical interpretation of these ^ . 

The starting point for both ekpyrotic and cyclic cosmology is the brane picture 
of our universe as suggested by string theory. The various string theories and 
11-dimensional supergravity are related to each other by a web of dualities 
[108], and so any particular one of these theories could be used as a model 
for the early universe. However, certain theories have proved more intuitive 
than others in tackling particular problems, and it is the duality between 
10-dimensional heterotic string theory with gauge group Eg x £^8 and 11- 
dimensional supergravity [43,44] that has inspired the early universe models 
described in this review. Neglecting fermions, the action of the maximal, 11- 
dimensional supergravity is given by [31] 



jy 1 A AMNPQ 

Ml. ..Mil A r r 

"l728 ^MiM2M3^M4...M7^M8...Mii 



(201) 



^ Theorems about past-incompleteness [15] as well as entropy bounds [37] do not 
apply, because of the singular crunch/bang transitions [96]. 

^ Many early studies dealing with the embedding of ekpyrotic/cyclic models in 
string theory/supergravity have considered the space in between the branes to be 
a slice of Anti de Sitter spacetime, see e.g. [48,93,99,78]. Obtaining 5-dimensional 
AdS in supergravity is straightforward [88], but as soon as the branes are added, 
the story becomes much more complicated, and not all issues have been resolved at 
present [25,33,32,63]. This is why we only present the theoretically better-motivated 
embedding into heterotic M-theory here. 
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where we denote 11- dimensional quantities with hats. Kn determines the 
strength of gravity in 11 dimensions. Here Gmnpq = 4!(9[m^jvpq] is the 4- 
form field strength associated with the 3-form Amnp- The /i(3)G'(4)G'(4) term 
is a Chern-Simons term and it is topological in nature. It will be unimportant 
in what follows. 11-dimensional supergravity is a classical theory, which no- 
body knows how to quantize. However, we can calculate anomalies (which do 
not depend on the microscopic structure of spacetime), and these allow us a 
glimpse of quantum M-theory. If we consider the 11-dimensional spacetime to 
take the form of a 10-dimensional spacetime plus a line segment, then anomaly 
cancelation requires us to add the following boundary terms [43,44] 



^boundary g^^2 
1 
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(202) 



The gauge fields correspond to the group Eg and the tri?^ terms are required by 
super symmetry. In the limit that the line segment becomes vanishingly small, 
the two boundaries sit on top of one another, and we recover the supergravity 
approximation to the 10-dimensional Eg x Eg heterotic string theory. In fact, 
the size of the fine segment is related to the couphng constant of the heterotic 
theory (roughly because the size of the line segment is reinterpreted as the 
dilaton field in 10 dimensions), via the famous relation [108]: 



2 



Ru = gi (203) 

where Ru is the radius of the 11th dimension and gs is the string coupling 
constant. Thus, at strong couphng, the world looks 11-dimensional. As the 
boundary branes collide, the 11th dimension disappears, but at the same time 
the coupling goes to zero. Hence one might expect the collision to be rather 
mild. 

In order to get a realistic model of particle physics, we can now consider 
curling up 6 of the spatial dimensions into a Calabi-Yau manifold. Then it 
turns out that in order to obtain the right magnitude of Newton's constant in 4 
dimensions, the line segment must be larger than the radius of the Calabi-Yau 
manifold by a factor of about 30 or so [8,109]. Hence, going up in energy, the 
world looks 4-dimensional, then 5-dimensional (as the line segment becomes 
visible), and only above the GUT scale 11-dimensional. This is the reason for 
choosing a 5-dimensional braneworld setup in describing the dynamics of the 
universe close to the big bang. The relevant 5-dimensional theory is called 
heterotic M-theory, and it can be obtained by dimensionally reducing the 11- 
dimensional theory described above on a 6-dimensional Calabi-Yau manifold. 
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as follows [72,73]. The metric is reduced according to 



d^n = VcY^'grnn^x^'^x'' + VcV^^bdxMx" . (204) 

where Vcy denotes the volume of the Calabi-Yau, gab being its metric. Indices 
m,n, ... run over 5 spacetime dimensions (including the line segment direction, 
labeled 11, but where we also use the coordinate x^^ = y), and a,b, ... denote 
Calabi-Yau directions. Due to the boundary terms (202), there is a modified 
Bianchi identity for the 4-form mode, which has to satisfy 



(dG)n««P« = {jty'{j''>Hy - 1) + J'^y + (205) 

where the currents are given by 

J{^) = tr^(i) /\ ^« _ 1 1-^^ /\ Jl (206) 

and where wc have chosen the coordinate positions y = ±1 for the location 
of the boundary branes. In the weakly coupled heterotic string theory the 
Bianchi identity reads 

dC ~ trR A R - trF(^) A F^^^ - trP^^^ A F^^^ . (207) 

The so-called "standard embedding" of the gauge connection into the spin 
connection ensures that 

trF(^) A F(^) = trR a R (208) 

and together with setting F^^^ = we see that we can have dG — and thus 
set G to zero consistently. In our case, there appears to be no way of achieving 
trp(i) /\ p(i) = itrR A R, and thus there will always be some components of 

G that are required to be non-zero. This mode is called the non-zero mode of 
the reduction. The standard embedding (208) then leads to 

{dG)nMNPQ = [^) {S{y-l)-S{y + l)){trRAR}^,^pQ. (209) 

The interpretation of this equation is that the ±tr7? A R terms are sources 
for oppositely charged branes sitting at the endpoints (y = ±1) of the line 
segment. We can solve the modified Bianchi identity (209) above (and at the 
same time satisfy the equation of motion for G) to lowest order by setting ^ 

Gabcd = ^^^bcd^ ^ef 0{y) , (210) 



2/3 



^ Compared to [72,73], we have re-scaled a such that a = aLOSw/3\/2- Note that 
the parameters a and /? used in this section are unrelated to the use of the same 
letters as indices of Hankel functions in section 4. 
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where 



e{y) = C^ for-l<y<+l mod 4 

^ -1 for - 3 < y < -1 mod 4 ^ ' 



7r fKii\l 



and 

where P = — Jc tri? A i? is the first Pontryagin class of the Calabi-Yau (the 
integration being done over the (2,2)-cycle C corresponding to the Kahler form 
u!ab)- Note that Gabcd is related to the integer /? and thus G is quantized; hence 
a is also quantized. We are only considering the lowest mode in the (Fourier) 
expansion of the full solution to the modified Bianchi identity, as all higher 
order modes are massive (their mass being of the order of the Calabi-Yau 
volume, i.e. the GUT scale) and decouple in the approximation considered 
here. In contrast to the other modes arising in the dimensional reduction of 
G, the non-zero mode cannot be set to zero. 

Keeping only gravity, the non-zero mode and the Calabi-Yau volume modulus 
in 5 dimensions (which is a consistent truncation), we end up with the action 



24 



5d 



1. 

2 



R - ^Vc^dmVcYd^VcY - 6a^Vc^ 



2^5 I J4d,y=+1 J4d,y=-1 ) 



(213) 



One should especially note the presence of the potential term Ga^/V^y which 
arose from the non-zero mode. This potential means that (l,4)-dimensional 
Minkowski space is not a solution to our theory, and hence not the vacuum. 
However such a "cosmological" term is dual in 5 dimensions to a 5-form field 
strength, which can support a 3-brane solution. In 5 dimensions a 3-brane is 
of CO dimension 1 and thus represents a domain wall. Note that in general for 
a 3-brane solution to be completely consistent, one needs a (4-dimensional) 
source action. But we actually have two 4-dimensional sources in the action 
(213) located at y = ±1. This double domain wall vacuum was found in [72]. 
It is given by 



ds^ = h''/'{y) [A^ i-dr' + df 2) + dy^] , 
VcY^B0'{y), 

h{y)^5ay + C, (214) 

where A, B and C are arbitrary constants. The y coordinate is taken to span 
the orbifold (line segment) 3^/2.2 with fixed points at y = ±1. In an 'extended' 
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picture of the solution, obtained by Z2-reflecting the solution across the branes, 
there is a downward-pointing kink aX y = —1 and an upward-pointing kink 
at 7/ = +1. These ensure the Israel conditions are satisfied, with the negative- 
tension brane being located at y = —1 and the positive-tension brane at 
y — +1. In the static domain wall solution above, the volume of the Calabi- 
Yau manifold and the distance between the boundary branes are determined 
in terms of the moduli B and C, while the scale factors on the branes are 
determined in terms of A and C. The modulus C additionally determines the 
height of the harmonic function /i at a given position in y. Before continuing by 
describing how we can extend this solution to allow for motion of the branes, 
we should add some remarks: by considering only the boundary branes, we are 
looking at a universal, but simplified version of heterotic M-theory. Typically, 
there are many more fields that can play a role, and there can be additional 
branes in the bulk. Such extra ingredients are usually necessary in order to 
obtain a realistic model of particle physics on the branes, see e.g. [16,21]. 
Here wc focus on the simplest case of having only boundary branes, since we 
are interested in the collision of these boundary branes. Also, the braneworld 
picture suggest a natural place for dark matter to reside, namely on the second 
boundary brane. Any form of matter residing on the "other" brane will appear 
as dark to us, since we would only detect its effects via gravity. However, there 
is no quantitative argument known for how much of the dark matter should be 
ordinary matter on second brane, or weakly-interacting matter on the visible 
brane. 

We are interested in finding the effective 4-dimensional theory describing the 

motion of the boundary branes that we have just discussed. We will obtain this 
by the method of the moduli space approximation [60]: we promote the moduli 
A,B,C in the solution above to arbitrary functions of the brane conformal 
time T, yielding the ansatz: 



Let us give a brief justification for this ansatz: firstly, we note that the ansatz 
satisfies the ry Einstein equation identically. This is important, since other- 
wise the ry equation would act as a constraint, see e.g. [40]. Secondly, there is 
no gT-y modulus, since this metric component is odd under the Z2 symmetry, 
and therefore has to vanish at the location of the branes. Any such component 
which is zero at the location of the branes, but non-zero in the bulk, is nec- 
essarily massive. In fact, from the work of [62], we know that, apart from the 
above moduli, all other perturbations have a positive mass squared. Having 
defined the time-dependent moduli, we would now like to derive the action 
summarizing their equations of motion. This is achieved by simply plugging 



ds' = h'^'{r,y) [A\r){ 
Vcy = B{t) h'l\T,y), 



dr^ + df2) + 52(r)d|/2 , 





(215) 
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the ansatz (215) into the action (213), yielding the result 



Smod = -6 / A^BIs [ 



A) 12\B J AB 

^-C" + -^—-i (216) 
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where we have defined 



4= /' /i" = --i--[(C + 5a)("+i) -{C- (217) 
j-i 5Q;(n + Ij 

This action can be greatly simplified by introducing the field redefinitions 



a^ = A^BIi, 

5 



02 
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(218) 
(219) 

(220) 



Note that a has the interpretation of being roughly the four-dimensional scale 
factor, whereas 0i and 02 are four-dimensional scalars. The definition (220) 
can be rewritten as stating that 



d02 = 



Vedc 



2(C + 5q;)V5 (C-5q;)V5 /a' 

5 

This expression can be integrated to yield 

'(1 + g2v/273</.2-)5/4 + (1 _ g2y^fli2-)5/4 
(1 + g2-^/273</'2)5/4 - (1 - e2"\/273<^2)5/4 



(221) 



C = 5a 



(222) 



In terms of a, 0i and 02, the moduli space action (216) then reduces to the 
remarkably simple form [60] 



^mod= / [-6a'Va^(0? + 0;^)]. 

J4d 



(223) 



The minus sign in front of the kinetic term for a is characteristic of gravity, 
and in fact this is the action for gravity with scale factor a and two minimally 
coupled scalar fields that is used for the calculations of cosmological pertur- 
bations in section 4. There is also a manifest 0{2) rotation symmetry for the 
scalar fields. The equation of motion for a reads 



6a" = -a (0? + 0^2), 



k'2\ 



(224) 
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while the equations of motion for (f)i and 02 immediately lead to the conserved 
charges Qi and Q2, according to 

aVi = <3i, aV2 = <92. (225) 
The solutions to these equations are given by 



a' = 2^Ql + Ql (r-r„), (226) 



iQi 



X- ln[02,o(r-r„)], (228) 



where Qi, Q2, Ta, 0i,o, and 02,o are constants of integration. 

We can now return to the ansatz (215) and relate physical quantities in five 
dimensions to the moduli fields a, 4>i and if we denote the distance between 
the branes by d, and the volume of the Calabi-Yau and the brane scale factors 
at the locations y = ±1 by Vcy± and a± respectively, then we have the 
relations 



d^ 



CY± ■ 



3(2a)V4 

(2c,)3/4g.Ai/v^ < 



g?ii/V2-072</)2 g2y^,^2)3/2 _ |l _ g^V^'^^ |3/2j^229) 

3/2 



(cosh ^27302^ 



1/4 



a, = (2a)VSae-^^/^^ 

(-sinhpl?>(t)2) . 



(230) 



(231) 



These relations are useful in interpreting particular solutions to the moduli 
equations of motion. Note in particular that the 4-dimensional effective theory 
describes the motion of both branes at the same time, as well as the geometry 
of the 5-dimensional bulk. In the limit that 02 — ^ +00, we have 



rf~(2a)-^/^e^^/^-^^/^. 



(232) 
(233) 



whereas for 02 — —00, we have 
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(234) 
(235) 



Thus, in both hmits, \nd and InVcyi are orthogonal variables. This means 
that, sufficiently far away from the 02 = axis, the fields 0i and 4>2 are, up 
to a re-scaling, simply related to In d and In Vcy± by a rotation in field space. 
Since the moduli space trajectories in terms of 0i and 02 are straight lines, 
far from the 02 = axis, the trajectories will also be approximately straight 
lines in terms of In d and In Vcy±- 

The 02 = axis plays a special role, in that it represents a boundary to moduli 
space. Positive values of 02 would lead to an unphysical negative volume of 
the Calabi-Yau manifold at the location of the negative-tension brane, and are 
forbidden. Physically, the boundary corresponds to a naked singularity in the 
bulk spacetime [67]. In fact, it can be shown that in the presence of a wide 
variety of matter fields on the negative-tension brane, the negative-tension 
brane bounces off this singularity in a smooth way and without encounter- 
ing the singularity [67]. This is only possible because of the special proper- 
ties of gravity on an object of negative tension. In comparing the moduli to 
higher- dimensional quantities via equations (229)-(231), it is apparent that 
this bounce of the negative-tension brane is equivalent to fiipping the sign of 
02 and thus 02 also, and hence the otherwise straight trajectories in scalar 
field space get reflected off the 02 = boundary, as shown in Fig. 6. This 
reflection has the important consequence of converting entropy perturbations 
around the background evolution into curvature perturbation, as discussed in 
section 4.1.2. 

As we have just seen, the general solution to the 4-dimensional effective theory 
(223) is a straight line of arbitrary slope in scalar field space, but bending 
at the boundary 02 = 0. In general, from the higher-dimensional point of 
view, the solutions correspond to a Kasner spacetime near the brane collision. 
However, one particular slope is special in that it corresponds to the least 
singular colliding branes solution studied in [59]. It is given by 



a = |2yor^/^ e^^ = |2yor|-^/^^, e^^ = |4a|/or|^/2^, (236) 



where 2yo has the interpretation of being the velocity of the branes at the col- 
lision, as, for small r, we have from (229) that d ~ 2yQr. For this solution, the 
brane scale factors and the Calabi-Yau volume are finite at the brane collision 
(as can be verified directly using (229)-(231)), and hence the spacetime metric 
asymptotes to Milne near the collision. It was argued in section 3 that we can 
take this solution as a model spacetime for the crunch/bang transition. As 
discussed there, the ekpyrotic phase ensures that the curvatures remain small 
until the quantum gravity regime is reached. In particular, higher-derivative 
corrections involving the Riemann tensor will also be small near the collision 
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Fig. 12. A Kruskal-type spacetime diagram in which the exact trajectories of the 
positive- and negative-tension branes are plotted. The negative-tension brane (in 
green) is always to the right of the positive-tension brane (in red), thereby shielding 
it from the naked singularity (thick black line). The bounces of the negative-tension 
brane off the naked singularity, as well as the collision of the branes themselves, 
are shown at a magnified scale in the inset. Note that the physical spacetime is 
restricted to the region bounded by the two branes. 

for the quasi-Milne spacetime, and hence we should be able to trust the above 
solution up to about a Planck time (or so) on either side of the collision. 

Note that the solution above is simply a solution to the moduli space effective 
action, and hence is only the leading order solution in an expansion in the 
brane velocity. However, the corresponding exact solution may actually be 
found directly in 5-dimensional heterotic M-theory (213) [59]. Above, we have 
used a gauge in which the bulk is dynamical and the branes are kept at fixed 
coordinate positions. Here, it will be more useful to use a gauge where we 
have a static bulk through which the branes are moving. Then, subject to the 
boundary condition that the spacetime asymptotes to Milne at the collision, 
one can derive a Birkhoff-type theorem which fixes the bulk metric to be given 
by [59] 

ds2 = -/dT2 + ^dr2 + r2dx2, 

/(r) = aV-/i/, VcY = r\ (237) 

where < r < ^JJi/ a, and the coordinate T is unrelated to the time coordinate 
appearing in the section 3.2. Physically, this solution describes a timelike naked 
singularity located at r = 0. The coordinates above do not cover the whole 
spacetime manifold, but the maximal extension may easily be constructed 
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following the usual Kruskal procedure [59]. The branes then move in this 
background spacetime according to the Israel junction conditions, and this 
motion is shown in Fig. 12 (for a detailed derivation see [59]). One finds that 
the induced brane metrics are indeed cosmological, with scale factors given by 

a±^{l±4y/JlT±y/\ (238) 

where t± is the brane conformal time, and we have re-scaled a± to unity at 
the collision, which is taken to occur at t± = 0. Upon setting = 0, we imme- 
diately recover the static domain wall solution (214), after a suitable change 
of coordinates. More generally, we require > to avoid the appearance of 
imaginary scale factors, is directly related to the brane velocity yo, and one 
can verify that the scaling solution (236) is recovered as an approximation at 
low velocity [60]. For > 0, starting from the brane collision, the two branes 
proceed to separate. However, while the positive-tension brane travels out to 
large radii unchecked, the negative-tension brane reaches the naked singular- 
ity (at which a_ and Vcy- tend to zero) in a finite brane conformal time 
T_ = (4y7^)~^. The resulting singularity is simple to regularize. As stated 
above, if almost any type of well-behaved matter is present on the negative- 
tension brane - even in only vanishing quantities - then, rather than hitting 
the singularity, the negative-tension brane will instead undergo a bounce at 
some small finite value of the scale factor and move away from the singularity. 
We can make this a little bit more precise by looking at the Friedmann equa- 
tion on the negative-tension brane. For the case in which a time-dependent 
scalar field (with a Vcy coupling to the Calabi-Yau volume scalar) is present 
on the brane, this equation takes the form (see [59]) 



where the constant p_ parameterizes the scalar kinetic energy density. The key 
feature of this equation is the negative sign in front of the first term on the 
right-hand side, refiecting the fact that matter on the negative-tension brane 
couples to gravity with the wrong sign. For sufficiently large values of the 
scale factor the right-hand side is dominated by the /xal^° term. If we further 
assume that the matter density on the branes is small compared to the brane 
tension (as is in any case necessary for the existence of a four-dimensional 
effective description, i.e. p_ <^ a), so that the term linear in p_ dominates 
over the quadratic term, then it follows that at some small value of the scale 
factor the entire right-hand side must vanish. Thus a negative-tension brane, 
initially traveling towards the singularity, will generically undergo a smooth 
bounce at some small value of the scale factor. After this bounce the brane 
travels away from the singularity back towards large values of the scale factor. 
This behavior is specific to the negative-tension brane, and moreover, persists 
even in the limit in which p_ (and hence the initial matter density) is negligibly 
small. It can also be generalized to a wide class of negative-tension branes and 
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matter types [67] . 



In order for the cyclic model and the ekpyrotic mechanism to work, we must 
assume a force between the boundary branes expressed by a potential of the 
form shown in Fig. 10. Could a potential of this shape arise from the heterotic 
M-theory setting under discussion here? First of all, wc should note that a neg- 
ative potential is completely natural in supcrgravity (hence all the problems 
in trying to explain the observed dark energy). Further, one can imagine that 
an attractive force between the branes could arise from virtual exchange of 
membranes stretching between the boundary branes (see [66] for an example of 
such membrane solutions). Preliminary studies do indicate that an attractive 
force can arise from open membrane instantons [80,68] in a certain parame- 
ter range. Unfortunately, a complete calculation of the potential involving all 
non-perturbative effects has been impossible to do so far. The vanishing of the 
potential near the brane collision seems natural, as the string coupling con- 
stant Qs goes to zero there. Hence one would expect non-perturbative effects 
to vanish also as ~ e~^/^^ [93]. Certainly, one would eventually like to sec a 
complete derivation of the potential in heterotic M-theory or a related setup. 



7 Conclusions and Outlook 

Postulating a slowly contracting phase of the universe with equation of state 
w ^ 1 (preceding the standard expanding phase of big bang cosmology) has 
extremely powerful consequences: it renders the universe flat and homogeneous 
over large patches, it delays the onset of chaos until the quantum gravity 
regime is reached, it prepares the universe ideally for a big crunch/big bang 
transition, and, on top of that, it generates small-scale scalar perturbations 
with a nearly scale-invariant spectrum. Hence, even if the higher-dimensional 
picture and the link to string theory seem unappealing to you, at the very 
least the ekpyrotic phase represents a new mechanism for solving the standard 
cosmological puzzles and a new way of generating cosmological perturbations. 
Moreover, these perturbations contain a significant non-gaussian component, 
so that they can be tested observationally and compared to simple, gaussian, 
inflationary models. Further into the future, the prediction of a blue small- 
amplitude gravitational spectrum will also be tested against observations. 

But more generally, ekpyrotic and cyclic models present a very broad con- 
ceptual framework, in which the big bang is not the beginning of time, but 
rather a physical event amenable to a testable physical description. Cyclic 
models allow for evolutionary timescales much larger than the time since the 
last big bang, and can provide qualitatively new explanations in this way. The 
embedding of ekpyrotic/cyclic models in string theory provides a geometric 
interpretation for many cosmological processes, and can eventually provide 
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observational tests of the theory. Many open issues remain, which should be 
tractable within string theory: in particular, the issue of moduli stabilization 
requires further work, as does the closely related question of "initial" condi- 
tions in the presence of many moduli. Also, the all-important cyclic potential 
has not yet been derived from first principles. Most promising is the recent 
realization that the brane collision, and the corresponding 4-dimensional big 
crunch/big bang transition, can perhaps be understood within the context of 
generalizations of the non-perturbative gauge/gravity correspondence. Much 
work remains to be done here before a physically realistic case can be treated; 
but this is a topic of such importance to cosmology that we should see this 
as an opportunity, rather than a drawback. It would be truly magnificent if 
we could find observational signatures of higher-dimensional physics, and of a 
time before the big bang. 
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